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•\  KJ'&L  Q*>jV*  •  ABSTRACT 

Topological  methods  may  be  employed  in  an  effective  manner  in  the 
analysis  and  design  of  a  communication  net.   One  application  is  in  the 
determination  of  the  reliability  of  the  net  with  respect  to  both 
complete  communication  and  k-terminal  communication.   The  effect  on  this 
reliability  caused  by  the  removal  of  one  link  is  studied. 

The  k-terminal  reliability  function  is  defined  and  its  application 
to  analysis  and  design  of  communication  nets  is  demonstrated.   A  digital 
computer  program  is  presented,  and  examples  of  its  use  are  included. 
A  formula  is  given  to  find  the  value  of  the  2 -terminal  function. 
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1.    Introduction. 

Any  system  that  provides  for  communication,  be  it  the  transmission 
of  information  or  the  transportation  of  commodities,  is  subject  to  pos- 
sible failure.   The  loss  of  security,  the  malfunction  of  a  radio  trans- 
mitter, and  the  destruction  of  a  bridge  are  examples  of  this  failure. 
Linear  graph  theory,  as  described  in  texts  such  as  [5],  has  been  ap- 
plied to  measure  the  reliability  of  the  system  given  some  measure  of  the 
reliability  of  the  parts.   Fu  and  Yau,  [4],  have  indicated  a  value  for 
the  probability  that  two  stations  can  communicate;  and  a  digital  computer 
algorithm  for  this  has  been  given,  [6].   The  probability  that  all  pairs 
of  stations  can  communicate  within  a  given  time  interval  has  been  studied 
by  Chan,  [2],  and  Fu,  [3]. 

In  this  work  the  application  of  the  k-terminal  reliability  function, 
to  be  defined  in  Section  3,  is  presented,  and  a  digital  computer  pro- 
gram is  included.   In  addition,  the  application  of  the  power  of  this 
function  without  computer  assistance  is  demonstrated.   A  formula  is  given 
to  find  the  value  of  the  2-terminal  reliability  function  without  enumer- 
ating the  Table  of  Combinations. 


2.   Statement  of  the  Problem. 

A  communication  net  may  be  represented  by  a  linear  graph  in  which 
the  edges  of  the  graph  represent  the  communication  links,  and  the  nodes 
of  the  graph  represent  the  stations  of  the  communication  net. 

A  weight  can  be  assigned  to  each  edge  of  the  graph  to  indicate  the 
probability  that  the  particular  communication  link  represented  by  the 
edge  will  operate  successfully.   This  weight  is  called  the  reliability 
of  that  edge.   This  probability  is  not  a  conditional  probability.   Each 
edge  is  assumed  to  be  associated  with  an  independent  random  variable, 
having  only  two  possible  states:   normal  operation,  and  not  operational. 

The  communication  net  can  be  analyzed  by  studying  the  topology  of 
the  graph.   Topological  methods  can  be  applied  on  several  different 
problems.   For  example,  if  one  link  of  a  communication  net  has  been  lost 
(destroyed),  it  must  be  decided  whether  or  not  to  replace  this  link  with 
one  of  the  other  existing  links  in  the  net.   The  information  needed  to 
make  this  decision  must  include  the  relative  importance  of  each  link, 
some  measure  of  the  quality  of  the  branch  (e.g.,  its  reliability),  and 
the  role  of  the  branch  in  the  net.   The  analysis  of  the  data  must  in- 
clude the  degree  of  importance  attached  to  the  overall  net  reliability 
as  contrasted  to  the  reliability  of  communication  between  a  pair  or  a 
group  of  terminals. 

Several  investigators,  T3],  [4],  [6],  have  proposed  quantitative 
measures  for  the  reliability  of  a  communication  net.   Investigations  by 
Chan  [2],  have  led  to  a  result  which  uses  the  reliability  function  as  a 
k-terminal  reliability  function.   This  function  may  be  applied  to  the 
reliability  of  communication  between  any  k-number  of  nodes.   The  topo- 
logical methods  which  are  used  to  generate  the  function  are  suitable 


for  use  on  a  digital  computer.  The  large  storage  requirement  is  the 
limiting  factor  in  their  application. 


3.   Properties  of  Reliability  Functions. 
3.1  2-terminal  Reliability  Function. 

The  reliability  function  is  defined  in  reference  [2].   Some  of 

the  most  important  definitions  are  included  here  for  ready  reference  and 

emphasis.   Using  the  notation  of  [2],  there  is  a  value  of  reliability, 

p.,  p0,...p  ,  assigned  to  each  of  the  edges  of  the  graph,  e  ,  e  ,...e  , 
1   2     e  i   2     K 

...e  .   Each  branch  of  the  net  is  assumed  to  be  associated  with  an  in- 
e 

dependent  random  variable.   If  the  branch,  b,  ,  is  operational,  the  edge, 

e.,  of  the  graph  assumes  the  value  p^ ,  which  is  the  probability  that 
i  i 

edge,  e.,  will  operate  normally.   If  the  branch  is  not  operational,  the 
edge,  e  ,  assumes  the  value  p\ ,  which  is  the  probability  that  the  edge, 

e  ,  will  not  be  operational.   Although  these  reliabilities  are  probabili- 

i 

ties,  they  are  considered  not  to  be  conditional  probabilities.   That  is, 
we  have  the  following  relations: 

Pi  +~±   -  1  (3.D 

Pr  (A  and  B)  -  Pr(A)Pr(B)  (3.2) 

Definition  3.1:   2-terminal  Reliability  Function: 

A  reliability  function  expressed  in  terms  of  the  edge  reliabilities  of 
a  network  is  a  2-terminal  reliability  function  if  it  represents  the 
reliability  between  two  distinct  terminals  (i.e.,  the  probability  of 
communication  between  these  two  terminals) . 

Let  q    represent  the  2-terminal  reliability  function  between 

x>y 

terminals  x  and  £.   Then  we  have 


'x,y  "VlV^e  +lYVy--  Pe-lPe  +  "•  +  m2eplP2' '  *Pe       (3'3) 
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where  m. ,  m«, . . ,m  , . . .m  e  are  either  1  or  0  depending  on  whether  the  com- 
bination of  states  of  the  edges  provide  a  communication  path  between  ter- 
minals x  and  2«   If  tne  path  exists,  then  m  =  1;  if  not,  m  =0.  When 
the  function  is  expressed  as  in  equation  (3.3),  it  is  said  to  be  in  its 
canonical  form. 

Definition  3.2:   Path  Product. 

A  path  set  is  a  set  of  branches  of  a  communication  path  between  a 
pair  of  terminals.   A  path  product  is  the  product  of  the  variables 
associated  with  the  elements  in  a  path  set. 

Definition  3.3:   Path  Product  Term  or  Primary  Path  Product  Term. 
A  canonical  term  is  a  path  product  term  if 

(1)  its  unbarred  variables  are  those  of  a  path  product; and 

(2)  all  other  variables  in  the  term  are  barred  variables. 

Definition  3.4:   Secondary  Path  Product  Term. 

A  canonical  term  is  a  secondary  path  product  term  if  its  unbarred 
variables  correspond  to  those  elements  forming  a  path  set  plus  one  or 
more  not  in  the  path  set. 

Definition  3.5:   Table  of  Combinations. 

The  collection  of  all  nonzero  canonical  terms  of  a  2-terminal 
reliability  function  is  called  the  Table  of  Combinations.   That  is, 

all  the  terms  in  which  m.  takes  on  the  value  1. 

l 

Definition  3.6:  k-edged  Path. 

A  k-edged  path  is  a  path  of  k  number  of  edges  or  branches  which 
constitute  a  communication  path  between  a  pair  of  terminals. 
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Theorem  3.1: 


Let  P.  be  an  A-edged  path  between  terminals  x  and  y_,  in  a  graph 


containing  e  number  of  edges.  The  Table  of  Combinations  for  the  2- 

nctioi 
(e-A) 


terminal  reliability  function,  q   ,  will  include  a  primary  path  product 


term  for  path  P.,  and  2V    -1  secondary  path  terms. 
A 

Proof: 

From  definition  (3.3)  and  (3.5),  the  primary  path  product  term 

(e-A) 

must  be  present.   There  are  2      possible  combinations  of  the  (e-A) 

variables  not  specified  by  the  primary  term.   However,  definition  (3.4) 

requires  the  number  of  unbarred  variables  to  be  at  least  (A+l) .  There- 

(e-A) 

fore,  the  maximum  possible  number  of  combinations  is  2     -1. 

If  a  graph  contains  e  number  of  edges  and  n  number  of  vertices,  then 
definition  (3.1)  indicates  that  there  will  be  2  number  of  terms  in  the 

2-terminal  reliability  function.   Theorem  (3.1)  indicates  that  one  path 

(e-A) 

of  A  number  of  edges  will  produce  2      terms  in  the  Table  of  Combina- 
tions.  This  is  the  same  as  saying  that,  from  definition  (3.1),  m.  will 

(e-A) 
be  nonzero  in  2      terms  of  the  2-terminal  reliability  function. 


Lemma  3.1: 

If  no  path  exists  between  two  terminals  x  and  y,  then  the  2-ter- 
minal reliability  function  will  have  no  nonzero  terms. 

Proof: 

By  definition  (3.1)  all  the  m  values  will  be  zero,  since  there  is 
no  path  between  terminals  x  and  y.   There  will  be  no  entries  in  the  Table 
of  Combinations. 
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Theorem  3.2: 

If  a  graph  contains  e  number  of  edges  and  n  number  of  vertices,  and 

if  an  A-edged  path  is  the  only  path  between  vertices  x  and  y,  then  there 

(e-A) 

will  be  2      number  of  terms  in  the  2-terminal  reliability  function, 

Sc.y' 

Proof; 

(e-A) 
By  theorem  (3.1)  there  will  be  2V     number  of  terms  for  the  A- 

edged  path.   Since  there  are  no  more  paths,  lemma  (3.1)  shows  there  will 

be  no  additional  terms.   Hence,  the  theorem. 

Theorem  3.3: 

Let  P  be  an  A-edged  path  and  P  a  B-edged  path,  each  between  ter- 
minals x  and  y_,  in  a  graph  containing  e   number  of  edges.   The  primary 
path  product  term  of  path  set  A  can  not  be  a  secondary  path  product  term 

for  path  set  B  in  q 

x,y 

Proof: 

The  unbarred  elements  in  the  secondary  path  product  terms  must  ap- 
pear as  barred  elements  in  the  primary  path  product  term  of  the  path  set 
A,  by  definition  (3.3).   Hence,  the  theorem. 

Theorem  3.4: 

Let  P.  and  P  be  path  sets  between  terminals  x  and  y_.   Then  at 
least  one  secondary  path  product  term  for  path  set  A  will  be  identical 
to  one  for  path  set  B. 

Proof: 

From  definition  (3.4),  the  collection  of  secondary  path  product  terms 
for  any  path  set  will  include  the  case  where  all  variables  are  unbarred. 
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Thi8  will  be  common  to  all  collections  of  secondary  path  product  terms 
Hence,  there  will  always  be  this  term  present. 


Example  3.1: 


Path  Set  A   PA  =[ei? 
Path  Set  B   PR  '^^ 


2-terminal  reliability  function,  q   , expressed  in  canonical  form: 

x»y 


q 
nx,y 

-         mlPlP2P3 
m2P1P2P3 

m3PlP2P3 

+ 
+ 
+ 

V)    in 

(8) 

(9) 

>) 

W2P3 

+ 

(2) 

(ID 

m5Pl?2?3 

+ 

(1) 

m6PlP2P3 

+ 

(3) 

(12) 

m7PlP2P3 

+ 

(4) 

(13) 

mQPiPoP<* 

(5) 

(6) 

(14) 

(1)  primary  path  product  term  for  path  set  A 

(2)  primary  path  product  term  for  path  set  B 

(3) (4) (5)  secondary  path  product  terms  for  path  set  A 
(6)  secondary  path  product  term  for  path  set  B 

observe  that  (5)  and  (6)  are  identical 
(7) (8) (9)  m  terms  are  0  because  terms  do  not  contain  a  path 

between  terminals  x  and  y_ 
(10)  this  term  will  always  be  zero 

(11) (12) (13)  a  tree  of  the  graph  is  made  up  of  the  unbarred 
terms 
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(14)  a  tree  plus  a  chord  of  the  graph  is  made  up  by  the  un- 
barred terms 

3.2  Multi -terminal  Reliability  Function. 

Reliability  functions,  which  are  expressed  in  terms  of  the  edge 
reliabilities,  may  be  derived  in  order  to  represent  the  reliability  of 
transmission  between  more  than  two  terminals.   Thinking  of  a  communica- 
tion net,  one  can  envision  a  classification  of  the  terminals  into  the 
categories  of  transmitter,  receiver,  or  both,  and  repeaters.   The  re- 
peaters may  be  looked  upon  as  intermediate  steps  in  the  communication 
paths.   All  of  this  may  enter  into  an  evaluation  of  the  reliability  of 
simultaneous  communication  between  all  pairs  of  terminals. 

Definition  3.7:   Complete  Communication. 

The  event  that  every  pair  of  terminals  x  and  y_  in  a  net  can  communi- 
cate with  each  other  simultaneously  is  complete  communication;  the  simul- 
taneous communication  between  all  pairs  of  terminals  of  the  net. 

Definition  3.8:   N-terminal  Reliability. 

The  probability  of  complete  communication  in  a  net  with  N-terminals 

is  called  the  "N-terminal  reliability11  of  the  net  and  is  denoted  by  q  . 

^n 

Definition  3.9:   k-terminal  Reliability. 

In  a  set  of  n  number  of  terminals,  consider  k  number  of  terminals 
of  transmission,  where  k  <  n.  Then,  the  probability  of  complete  com- 
munication among  the  set  of  k  number  of  terminals  is  called  the  k-ter- 
minal reliability  of  the  net  and  is  denoted  by  q,  . 

k 

From  definitions  (3.7)  and  (3.9),  a  net  may  have  several  differ- 
ent values  of  k-terminal  reliability  depending  on  which  terminals  are 
included  in  the  set  of  k  number  of  terminals  of  transmission.   The 
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choice  of  terminals  included  in  this  set  is  of  fundamental  importance 
in  the  interpretation  of  the  analysis  or  data  collection. 

For  a  given  net,  q  is  found  in  the  following  way:   first  determine 
q    for  each  pair  of  vertices.   There  will  be  %(n)(n-l)  for  a  net  with 

x»y 

n  number  of  terminals.   Then  sum  all  the  canonical  terms  which  are  common 

to  all  q    functions.   To  find  q,  ,  find  q    for  each  pair  of  terminals 
Mx,y  ^k       nx,y 

in  the  set,  K.  There  will  be  %(k)(k-l)  number.   Then  sum  all  the  canoni- 
cal terms  common  to  all  these  q    functions. 

x»y 

Another  method  is  first  to  determine  a  complete  set  of  trees  for  the 
net.   Complete  communication  is  possible  if  and  only  if  the  unbarred  vari- 
ables of  a  canonical  term  include  those  corresponding  to  the  branches  of 
a  tree.  By  taking  the  sum  of  all  canonical  terms  that  contain  a  tree, 
or  a  tree  plus  chords,  made  up  from  the  unbarred  variables,  q  is  deter- 
mined.  If  the  terms  which  are  summed  contain  a  subtree  of  the  graph,  or 
a  subtree  plus  chords,  containing  all  the  vertices  of  the  set  K,  the  q, 
is  determined.   This  method  can  be  applied  on  the  canonical  form  result- 
ing after  the  first  q    calculation.  This  can  reduce  the  labor  required 

x»y 

to  find  q  or  q  provided  all  the  trees  of  the  graph  are  known. 

Path  set 


S.: 


*  2 
Figure  3.1  repeated 

Part  A 


mlPlP2P3  + 
m2p1p2p3  + 

m3P1P2P3  +  (1) 

m4P1P2P3  +  (3)  (7) 


Path  set   PD  -  $exe3^ 

-hi 

=  he2} 


Path  set   P„  = 


Path  set 


Part  B 


Vz  =  mlPlP2p3  + 


m2PlP2P3  +  ^ 
m3PlP2P3  + 
m4plP2p3  +  *3^  (7^ 
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m5PlP2P3  + 


m6P1P2P3  +  (2)  (8) 
m7plp2p3  +  (4)  (9) 


m8PlP2P3 


(5)  (6)  (10) 


m5PlP2P3  + 

m6plp2P3  +  (4)  (8) 

m7P1P2P3  +  (2)  (9) 

m8PlP2P3   (5)  (6)  (10) 


(1)  Primary  path  product  term  for  path  set  C  and  E. 

(2)  Primary  path  product  term  for  path  set  D  and  F. 

(3) (4) (5)  Secondary  path  product  term  for  path  set  C  and  E. 

(6)  Secondary  path  product  term  for  path  set  D  and  F. 

(7)  (8)  (9)  (10)  make  up  a  tree,  or  a  tree  plus  a  chord,  of 

the  graph,  made  up  of  the  unbarred  terms. 


Example  3.3 


Number  of  trees,   t 


D 


-  3 


Figure  3.1  repeated 

The  3  trees  of  the  graph  are: 

Using  the  results  of  examples  (3.2)  and  (3.1),  we  have 

%   =  W2p3  + 

m6Plp2p3  + 

m7Plp2P3  + 

m8PlP2P3 
by  the  process  of  summing  the  canonical  terms  in  common.   From  example 

(3.1)  alone,  or  (3.2)  part  (a)  alone  or  part  (b)  alone,  the  same  results 
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could  have  been  obtained  using  the  method  of  trees. 

In  this  simple  example,  for  any  set  of  terminals  of  transmission 
less  than  the  given  set,  the  problem  reduces  to  one  of  the  examples, 
since  it  becomes  a  2-terminal  reliability  problem. 

3.2  Some  Characteristics  of  the  Path  Product  Terms. 

A  graph  contains  5  edges.  A  matrix  can  be  constructed  in  which  the 
elments  are  the  number  of  path  product  terms  that  contain  U  number  of  un- 
barred variables  compared  to  the  L  number  of  edges  in  the  path. 

Number  of  unbarred  variables 
12    3    4    5 


Path  Length 
4 
3 
2 

1 


0 
0 
0 

1 


The  number  of  path  product  terms 
This  array  is  recognized  as  the  start  of  Pascal's  Triangle.   See 
reference  [1],   The  nonzero  matrix  elements  are  the  so  called  binary 
coefficients,  defined  by  the  formula 


(3.4) 


3/  J  -'(n-j).' 

where,  in  this  case,  n  =  (e-A),  and  A  is  the  path  length,  and  j  =  0, 
1,2,...,  (e-A). 

If  all  the  branches  of  the  net  have  the  same  reliability,  p,  then 
the  collection  of  all  path  product  terms  resulting  from  an  A-edged  path, 
in  a  graph  with  e  edges,  can  be  written  as  follows: 
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(e"A)  PAP(e-A)    +    (elA) 


(A+l)-(e-A-l)         /e-A]    <A+2).(e-A-2)     + 


P  P  +1     2    'p 


/e-A\      (e-l)-(l)        ,/e-A\  e 
'(e-A-lj  P  P  Ae"AJP  (3'5) 


e-A 


«D»(wl>p('4"l) 


i  =  0 


where  K(i)  is  the  collection  of  the  binary  coefficients  defined  by  the 
formula 

e-A  i 

n. 


i  J  i!(n-i).' 

Examples  (3.1)  and  (3.2)  illustrate  the  generation  of  the  primary  and 
secondary  path  product  terras.   Theorem  (3.4)  points  out  that  there  will 
always  be  at  least  one  of  the  secondary  path  product  terms  common  to 
the  expansions  generated  from  two  or  more  path  sets  is  a  2-terminal  re- 
liability function.  Theorem  (3.3)  proves  that  the  primary  path  product 
term  is  never  a  secondary  path  product  term  for  another  path  in  the  2- 
terminal  reliability  function;  that  is,  it  is  not  one  of  the  duplicate 
terms.   The  exact  number  of  unique  terms  added  to  this  function  is  a 
problem  in  conditional  probability.   The  following  group  of  theorems  is 
fundamental  to  the  solution  of  this  problem. 

The  theorems  and  lemmae  in  this  section  are  subjected  to  the  follow- 
ing conditions:   a  graph  contains  e  number  of  edges  and  n  number  of  verti- 
ces.  Path  set,  P.     ,  is  a  1-edged  path  set  between  terminals  x  and  y_, 
made  up  of  branch  b,  alone. 
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Theorem  3.5; 

No  multi-edged  path  between  terminals  x  and  y  may  contain  branch  b. 

Proof: 

By  definition  (3.2)  a  path  set  is  a  set  of  branches  of  a  communica- 
tion path  between  a  pair  of  termunals.   No  additional  branches  are  in- 
cluded in  the  path  set.   Hence,  since  branch  b  is  a  path  by  itself,  any 
set  that  is  made  up  of  branch  b  and  one  or  more  additional  branches  can- 
not be  a  path  set  between  terminals  x  and  y_.   Definition  (3.6)  of  a  k- 
edged  path  allows  the  conclusion  that  no  multi-edged  path  between  ter- 
minals x  and  y_  contain  branch  b. 

Lemma  3.2; 

There  can  be  no  primary  path  product  term  in  which  branch  b  ap- 
pears as  an  unbarred  variable  in  the  2- terminal  reliability  function, 

q   »  except  path  set  P, 
x,y'  l>xy 

Proof; 

Lemma  (3.2)  follows  immediately  from  theorem  (3.5)  and  definition 
(3.3). 

Lemma  3.3; 

There  can  be  no  secondary  path  product  terms  generated  from  any 

path  between  x  and  y  in  which  branch  b  appears  as  an  unbarred  variable 

except  path  set  P, 

l,xy 

Proof: 

This  lemma  follows  directly  from  theorem  (3.5)  and  definition  (3.4). 
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Theorem  3.6; 

Let  path  set,  ?0        ,  be  an  A-edged  path  set  between  terminals  x 
Z,xy 

(e-A-1) 
and  y_.  This  path  will  generate  2       secondary  terms  that  are  dupli- 
cates of  those  generated  by  P. 

1  ,xy 

Proof: 

The  number  of  possible  unique  combinations  of  the  (e-A)  number  of 

barred  variables  has  decreased  by  a  factor  of  2  because  the  1-edged 

(e-1) 

path,  P    ,  has  previously  formed  2      combinations  of  the  same  vari- 
1  ,xy 

ables.   Only  those  in  which  the  variable  b  appears  as  a  barred  variable 

(e-A-1)  (e-A) 

remain  ungenerated.   That  number  is  2       .   There  were  2      possible, 

(e-A-1)  (e-A-1) 

but  now  only  2       ,  meaning  that  2       were  duplicates. 

Investigations  have  shown  that  the  number  of  duplicates  generated 

among  all  the  path  product  terms  are  arranged  in  correspondence  with  the 

binary  coefficients.  That  is,  if  there  are  eight  duplicated  generated, 

there  will  be  a  1-3-3-1  distribution,  based  on  the  number  of  unbarred 

variables  in  the  path.   If  there  are  D  number  of  duplicates,  then  the 

distribution  will  be  as  given  by 

(j)*  irrb)  > J  =  0-  l>2 D 

The  number  of  unique  path  product  terms  that  are  generated  by 

each  path  set,  P.,  in  the  2-terminal  reliability  function,  q   ,  can 
i  x,y 

be  found  by  the  following  method. 

a)  Write  the  path  matrix,  order  (p  x  e) ,  where  £  is  the 
number  of  path  sets  between  terminals  x  and  y,  and 
£  number  of  edges  in  the  graph. 

b)  Rearrange  the  rows  such  that  the  shortest  path  is  row  #1, 
the  next  shortest  is  row  #2,  and  so  forth,  until  the  last 
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row  contains  the  longest  path. 
c)   Count  the  number  of  columns  in  the  row  in  which  there  is 
no  entry  in  any  of  the  preceeding  rows.  For  example, 
this  number  will  be  (e-k)  for  the  first  row.   This  num- 


ber may  be  zero.   Call  this  number  n  . 


d) 


Calculate  2  i.   This  result  is  the  number  of  unique  path 
product  terms  generated  by  the  path  set  in  the  2-terminal 
reliability  function. 
The  total  number  of  terms  in  the  Table  of  Combinations  is  the  sum 
of  the  results  of  step  (d)  above.   It  is  seen  that  if  n  is  zero,  that 
one  term  is  still  added  to  the  Table  of  Combinations  by  the  path,  which 
agrees  with  Theorem  (3.3).   If  there  is  only  one  path,  then  there  will 
be  2      number  of  terms  which  agrees  with  Theorem  (3.2).   Further 
statements  which  tend  to  prove  the  validity  of  this  method  follow  the 
presentation  of  this  example. 


Example  3.4; 


Figure  3.2 
Path  Matrix: 

edges : 
paths:  i 

2 
3 


The  paths  between  vertices  A  and  B: 


Wsj  •  fal  -WvJ 
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edges: 

1 

paths: 

2 

0 

1 

1 

3 

1 

After  rearranging,  the  path  matrix  becomes: 

2  3  4  5 
0  10  0 
0  0  0  1 
10        10 

Counting  the  number  of  columns  in  each  row  in  which  there  is  no  entry 

in  that  column  in  any  of  the  preceeding  rows: 

paths :    2  4 

1  2 

3  0 


Calculating: 


2ni: 


paths:    2  16 

1  4 

3  1 

From  definitions  (3.1  -  5),  there  are  established  conditions  which 
a  collection  of  variables  must  meet  in  order  to  be  an  entry  in  the 

Table  of  Combinations.   Looking  at  example  (3.4),  one  can  see  that  path 

4 
(2)  will  generate  2  unique  terms.  After  they  have  been  generated,  path 

3 
(1)  generates  a  collection  of  2  number  of  possible  terms.   All  these  in 

which  edge  (3)  appears  as  an  unbarred  variable  will  be  duplicates  of 

those  generated  by  path  (1).   This  reduces  the  number  of  terms  by  a 

2 
factor  of  2,  yielding  the  result:  2  .   This  can  be  understood  from  the 

concept  of  combinations;  there  were  only  2  elements  to  be  combined  in 

2 
2  different  ways  vice  what  had  appeared  to  be  3  elements;  i.e.,  2  vice 

23. 
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The  same  argument  holds  for  path  (3).  While  there  would  appear 
to  be  2  elements  that  can  be  varied,  there  are  none  which  have  not  al- 
ready been  varied.   The  fact  that  the  primary  path  product  term  is  unique 
itself  is  given  in  theorem  (3.3).   It  is  from  this  logic  that  the  above 
stated  procedure  is  derived. 

The  knowledge  of  these  numbers  leads  to  a  method  to  determine  the 
value  of  the  2-terminal  reliability  function.   With  equation  (3.5)  in 
mind,  a  new  equation  can  be  written  for  each  path  set  using  the  numbers, 
n   found  by  the  above  described  procedure.   Let  P    be  a  k-edged  path 

i  x»y 

between  terminals  x  and  y  in  a  graph  with  e  number  of  edges.   For  this 
path,  let  n  ■  N.   Then  one  can  write 


'Vj-'-V  (?lV--V 


-(N-l) 
rrv     +  .  .  . 


N  \   (N-l). 

I  N-l 


where 


(3.6) 


is  the  value  of  all  the  path  product  terms  added  to  the 


Table  of  Combinations  by  the  path,  P  . 

Where  the  first  factor  is  the  product  of  the  k  number  of  branch 

reliabilities  associated  with  the  branches  in  the  path,  P   ;  the  second 

x,y' 

factor  is  the  product  of  the  (e-k-N)  number  of  complements  of  the  branch 
reliabilities  of  the  branches  which  have  previously  appeared  as  unbarred 
variables  in  other  path  product  terms;   the  third  factor  is  the  summa- 
tion of  all  the  path  product  terms  in  which  N  number  of  variables,  associ- 
ated with  the  branches  not  in  the  path,  P   ,  and  not  previously  speci- 
fy 

fied  as  being  unbarred  in  other  path  product  terms,  appear  as  barred  and 
unbarred  elements. 
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However,  the  quantity  inside  the  brackets  is  equal  to  one.  This 
statement  is  true  because  the  quantity  inside  the  brackets  is  an  applica- 
tion of  the  binomial  theorem  for  two  conditions:   the  branch  reliabilities 
being  all  equal;  and  the  branch  reliabilities  being  not  all  equal. 
Case  1)   Let  all  the  branch  reliabilities  be  equal  to  r. 

Then  the  quantity  inside  the  brackets  of  equation 
(3.6)  becomes 


(r  +  r) 


N 


Since  (r  +  r)  =  1,  then  (r  +  r)  =1. 
Case  2)   Let  the  reliability  of  branch,  b  ,  be  given  by  r. . 
Then  the  quantity  inside  the  brackets  of  equation 
(3.6)  may  be  rewritten: 

(rx  +  rx)(r2  +  r,,)  .  .  .  (rN  +  rN)  =  (1)(1)...(1)  =  1 
Therefore,  equation  (3.6)  can  be  rewritten 


(riY'-rk)(?LV--V  = 


In  the  event  all  branch  reliabilities  are  equal, 


(3.7) 


k.(e-k-N) 
r  r 


(3.8) 
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4.   Topological  Analysis  of  the  Net. 

4.1  Collection  of  All  Possible  Path  Sets. 

Fundamental  to  the  analysis  of  any  communication  net  is  the 
knowledge  of  all  possible  paths  between  the  terminals  x  and  y_.   Since  no 
path  may  be  ignored,  even  in  relatively  small,  simple  nets,  the  task  of 
finding  all  paths  is  a  laborous  task.   Linear  graph  theory  provides  an 
efficient  method  to  accomplish  this  task. 

Consider  a  graph  with  e_  number  of  edges  and  n  number  of  vertices, 
and  define  a  circuit  to  be  a  closed  edge  train.   Then,  the  theory  of 
topology  has  shown  there  to  be  (e  -  n  +  1)  ■  w  number  of  independnet  cir- 
cuits in  the  graph.   These  are  called  the  fundamental  circuit  sets,  and 
are  denoted  by  Bf.  Any  additional  circuits  that  appear  to  exist  in  the 
graph  can  be  shown  to  be  algebraically  derived  from  the  set  of  indepen- 
dent circuits. 

Next  there  must  be  built  a  collection  of  sets,  called  set  B,  which 
includes  the  null  set,  the  W  number  of  sets  of  B_,  and  the  collection  of 
all  possible  combinations  of  circuit  sets,  which  are  generated  by  the 
ringsum  operation  of  two  or  more  circuit  sets  in  Bf .   The  ringsum  opera- 
tion is  defined  as  modulo  2  addition  of  the  respective  column  entries  of 

w 
the  two  sets.   (See  example, (4. 2) .   The  set  B  will  contain  2  number  of 

sets. 

A  path  set,  P   ,  is  chosen,  and  expressed  as  a  row  matrix  of  e_ 

columns.   Set  B  is  represented  in  matrix  form  of  order  (w  x  e) .   The 

ringsum  operation  of  P    and  each  set  in  B  produces  a  new  set  which  is 

x»y 

written  in  a  collection  of  sets  called  set  BRP.  This  is  also  in  matrix 
form  of  order  (w  x  e).  Linear  graph  theory  shows  that  the  ringsum  of  a 
path  set  and  a  circuit  set  will  be  either  a  new  path  set  or  the  disjoint 
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union  of  a  path  set  and  a  circuit  set.  Now  each  entry  in  BRP  must  be 
examined  to  retain  only  the  pure  path  sets,  which  are  written  in  set  PFC, 
also  in  matrix  form  of  order  (p  x  e) ,  where  £  is  the  number  of  paths 
between  x  and  y_. 

One  process  of  examining  set  BRP  is  to  look  at  each  entry  first  to 
see  whether  the  original  path  set,  P   ,  remains.   If  it  does,  and  if 

x»y 

there  are  other  branches  present,  then  it  is  not  a  pure  path  set  and  must 
be  discarded.   If  it  does  not,  then  it  must  be  ascertained  whether  any 
of  the  circuits  of  set  B  are  included.   If  so,  then  it  is  not  a  pure  path 
set.   Those  that  remain  are  the  pure  path  sets,  and  more  importantly, 
all  the  path  sets  between  terminals  x  and  y_.   A  block  diagram  of  this 
procedure  is  given  in  figure  (4.2). 

The  ringsum  operation  is  built  into  the  computer  program  by  means 
of  three  FORTRAN  IV  instructions.   The  block  diagram  for  this  operation 
is  shown  in  figure  (4.3).   Example  (4.1)  shows  the  ringsum  operation  for 
the  four  possible  cases.  The  parameters  must  be  logical  elements;  that 
is,  either  True,  or  False,  with  a  1  bit  representing  True,  and  a  0  bit 
False. 
Example  4.1: 

Case  1    Case  2    Case  3    Case  4 
A  =        0        1        1        0 
B  =        1        0        1        0 

C  =        1        1        0        0 

The  step-by-step  computer  program  solution  proceeds  as  follows: 

after  step  1,  C  =        0  0  0  0 

after  step  2,  C  =        1  1  1  0 

after  step  3,   C  ■        1  1  0  0 
The  final  results  are  identical. 
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Figure  4.1 


Set  Bf: 


edge  number 
Circuit  number  1 


e2e3e4 


Circuit  (3)  =     \  e  c 
circuit  (1)  and  (2)  as  follows 

1    1 


0 


Choose  P 


A,B 


M 


Set  B: 


e  =  5 
m  =  4 
w  =  5-4  +  l  =  2 


2W  =  22  =  4 


Circuit  (1)  \ex  e2  e^i 
Circuit  (2)  £e3  e^   e5? 


3    4    5 
0    0    1 

Order  (w  x  e  =  2  x  5) 

is  derived  algebrically  by  ringsum  of 


edge  number 
set  number 

1 

1 

~b 

2 
0 

3 
0 

4 
0 

5 

"o 

2 

1 

1 

0 

0 

1 

3 

0 

0 

1 

1 

1 

4 

1 

1 

1 

1 

0 

Order  (2  x  e  =  4  x  5) 


PAB  ®  B  =  BRP 
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Set  BRP 


edge  number 
set  number 


1 
2 
3 

4 


2 

0 

1 
0 

1 


,w 


Order  (2  x  e  -  4  x  5) 


Set  PFC  = 


edge  number 
path  number 


Paths  are 


Order  (p  x  e  ■  3  x  5) 


K[  '  £e2  e5^  •   ^e2e3e4^ 


29 


calculate 
(e-n+1)  =  w 


write  w  inde- 
pendent circuits 
sets  in  B„ 


generate  set  B 
from  the  possi- 
ble combinations 
of  sets  in  B„ 


I 


choose  a  path 
set,  P 


i 


Ringsum  P 

and  B  to  obtain 
set  BRP 


I 


Find  pure  path 

sets  in  BRP; 
write  them  in 
set  PFC 


Block  Diagram  of  Process  to  Find  All  Paths 

Figure  4.2 
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sua. 


C  =  True 


B       ^ 

»vTru 

L   -   False 

alse 

\m 

Block  Diagram  of  Ringsum  Operation  in  Computer  Program 

Figure  4.3 
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4.2  Path  Finding  Program. 

A  digital  computer  program  can  be  devised  to  find  all  the  paths 
between  terminals  x  and  y_.   The  program  discussed  in  Appendix  (6)  ac- 
complishes this  task.   (For  further  details,  see  that  section). 

The  generation  of  all  possible  combinations  of  the  fundamental  cir- 
cuit sets  is  an  important  step,  and  no  one  may  be  omitted.   The  method 
used  by  the  program  to  find  them  all  is  the  following:   to  count  from  4 
to  7  in  binary,  one  adds  "100"  to  the  binary  numbers  from  0  through  3. 
Similarly  to  count  from  8  to  15,  one  adds  "1000"  to  the  binary  numbers 
from  0  to  7.   See  example  (4.2).   This  procedure  guarantees  that  all 
possible  combinations  are  obtained,  because  all  numbers  are  obtained  in 
the  analogy. 
Example  4.3; 


Decimal 

Binary 

0  = 

0 

0 

0 

0 

1  = 

0 

0 

0 

1 

2  - 

0 

0 

1 

0 

3  - 

0 

0 

1 

1 

to  count  from  4  to  7  in  binary,  add  100  to  those  numbers  preceeding  it. 


4  = 

0 

1 

0 

0 

5  - 

0 

1 

0 

1 

6  = 

0 

1 

1 

0 

7  = 

0 

1 

1 

1 

Similarly, 

from  8  to  15, 

add  1000 

8  = 

1 

0 

0 

0 

9  = 

1 

0 

0 

1 

10  = 

1 

0 

1 

0 

11  = 

1 

0 

1 

1 
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Decimal 

Binary 

12  ■ 

1 

1 

0 

0 

13  = 

1 

1 

0 

1 

14  = 

1 

1 

1 

0 

15  = 

1 

1 

1 

1 

4.3  Branch  Loading. 

When  a  communication  net  is  analyzed,  each  link  is  examined  in- 
dividually.  In  addition  to  the  reliability  of  the  link  to  perform  its 
function,  there  is  needed  a  measure  of  the  relative  importance  of  the 
link  in  the  various  modes  of  operation  of  the  complete  net.   For  example, 
the  importance  of  a  link  may  be  different,  if  complete  communication  is 
required,  than  if  only  K-terminal  communication  is  required.  One  such 
measure  of  relative  importance  of  a  link  is  called  the  load  of  a  branch. 

Fu  has  defined  the  load  of  a  branch  to  be  the  total  number  trees 
that  contain  this  branch,  [3].   As  shown  before,  complete  communication 
is  possible  only  when  the  net  contains  at  least  a  tree  in  the  operational 
state.   Since  each  branch  is  contained  in  one  or  more  trees  of  the  net, 
the  number  of  trees  that  contain  this  branch  can  be  found,  and  this  number 
used  to  indicate  the  relative  importance  of  the  branch  in  complete  com- 
munication.  This  number  can  be  normalized,  if  desired,  since  the  total 
load  of  all  branches  can  be  quickly  calculated  by  a  formula  shown  by  Fu, 
[3],  as  n 

X    h       =  <v  -  1>t 

i=l 

where  1.  (i  ■  1,  2,  . . . ,m)  is  the  load  of  branch  b. 
1  x 

in  a  network  with  v  number  of  vertices,  n  branches,  and  t   number 
of  trees. 
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This  method  of  branch  loading  determination  can  be  extended  to 
the  k-terminal  communication  case  by  means  of  subtrees.   However,  branch 
loading  can  also  be  defined  in  another  way. 

Definition  4.1:  Branch  Loading. 

The  load  of  a  branch,  b.,  is  defined  to  be  the  total  number  of  paths 

1 

(path  sets)  that  contain  this  branch. 

A  path  set  was  defined  to  be  a  set  of  branches  of  a  communication 
path  between  a  pair  of  terminals.  When  complete  communication  is  con- 
sidered, definition  (4.1)  makes  use  of  all  the  paths  between  all  pairs 
of  terminals.   This  means  that  the  numbers  obtained  for  the  various 
branch  loadings  will  be  larger  than  those  obtained  using  the  definition 
of  Fu;  however,  the  relative  magnitudes  will  be  the  same.  That  is,  the 
branch  with  the  highest  loading  will  be  identified  to  be  the  same  by 
either  definition.   The  advantage  of  definition  (4.1)  lies  in  the  simple 
way  to  evaluate  rather  than  having  to  find  the  trees  of  the  net. 

4.4  Path  Length 

When  a  communication  net  is  analyzed,  all  the  paths  between  ter- 
minal pairs  are  found.   The  loading  of  the  branches  is  measured.   Now 
the  paths  themselves  must  be  examined.   The  length  of  the  communication 
path  is  of  great  significance.  We  know  from  experience  that  the  more 
often  data  is  handled,  the  greater  the  chance  for  error.   This  is  con- 
firmed by  this  simple  example. 

Example  4.4: 

Let  p  be  the  reliability  of  all  branches,  b  ,  and  p<  1.   Then 

the  reliability  through  one  branch  is  p.  The  reliability  of  communica- 

2 
tion  through  two  branches  is  (  p  x  p)  =  p  ;  through  j  number  of  branches: 
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p  .   Since  p  <  1,  then  pJ <  p  for  all  j  >  1. 

This  effect  of  path  length  appears  directly  in  the  2-terminal  re- 
liability function  by  means  of  the  secondary  path  product  terras.   Section 
(4.6)  discusses  the  effect  of  the  secondary  path  product  terms.   See 
also  section  (5),  Application  of  Ideas  without  the  Aid  of  a  Computer. 

4.5  Branch  Reliability. 

Much  work  is  currently  being  done  by  many  investigators  in  the  fields 
of  Operations  Research,  Reliability,  Non -Destructive  Testing,  and  many 
others,  to  provide  a  high  quality  measurement  of  the  reliability  of  a 
communication  link.   This  work  does  not  concern  itself  with  the  finding 
of  this  value  of  reliability.   It  must  be  pointed  out,  however,  the  signi- 
ficance of  accuracy  in  this  value  and  its  effect  on  the  interpretation 
of  the  results. 

Results  have  shown  that  considering  a  net  with  e  =  10  branches,  and 
assuming  all  branch  reliabilities  to  be  equal  to  .90,  the  fourth  signi- 
ficant figure  has  risen  to  the  importance  level  of  being  the  deciding  fac- 
tor in  net  configuration. 

4.6  Secondary  Path  Product  Terms. 

The  secondary  path  product  terms  of  definition  (3.4)  are  not  of 
secondary  importance  in  the  determination  of  the  reliability  function. 
As  can  be  seen  from  example  (3.1),  the  majority  of  nonzero  canonical 
terms  are  secondary  path  product  terms.   The  exact  number  of  these  terms 

was  given  in  theorem  (3.1),  namely,  for  an  A-edged  path,  there  will  be 

(e-A) 
2     -1  number  of  secondary  path  product  terms. 

The  theorems  in  this  section  are  subjected  to  the  following  condi- 
tions: A  graph  contains  e  number  of  edges,  and  n  number  of  vertices. 
The  branch  reliability,  p   of  branch  b  ,  is  greater  than  its  complement, 
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p  ,  and  is  the  same  for  all  branches. 

Definition  4.2:  Value  of  a  path  product  term. 

The  value  of  a  path  product  term  is  defined  to  be  the  product  of 
the  branch  reliabilities  of  the  unbarred  variables  and  the  complement 
of  the  branch  reliabilities  of  the  barred  variables  of  the  path  pro- 
duct term.   It  is  applicable  to  both  primary  and  secondary  path  product 
terms.   It  always  takes  on  a  value  less  than  one,  since  it  is  a  product 
of  numbers  each  less  than  one. 

Theorem  4.1: 

Any  one  of  the  secondary  path  product  terms  has  a  larger  value  than 
the  primary  path  product  term  associated  with  it,  for  any  A-edged  path, 

Proof; 

The  value  of  a  primary  path  product  term  is  determined  using  de- 
finitions (3.3)  and  (4.2)  to  be 

A        -(e-A) 
P         P 

By  definition  (3.4)  a  secondary  path  product  term  must  have  at  least 
one  more  unbarred  variable  than  the  primary  term.   Therefore,  we  write 

(A  +  a)   -  (e-A-a) 

p       '  P 

where  a  is  a  positive  integer,  i.e.,  a  ^  1 

_-  A   -(e-A)    ^T      (A  +  a)     -(e-A-a) 

If  P  •  P       »  ^  P        .  Pv 

Then  1  <  pa  .   p("a) 

Since    P  <  p, 

-.       A      -(e-A) (A  +  a)      _(e-A-a) 

Thus    p    .   pv     <p      '.   pv      . 
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Corollary  4.1: 

The  sum  of  the  values  of  the  primary  and  all  the  secondary  path 
product  terms  of  a  1-edged  path  is  greater  than  that  of  any  multi-edged 
path  between  the  same  terminals  x  and  y_. 

Proof: 

This  Corollary  follows  directly  from  the  theorem  (4.1)  by  letting 
A  be  equal  to  1. 

Theorem  4.2: 

The  sum  of  the  values  of  all  the  primary  and  all  the  secondary  path 
product  terms  of  an  A-edged  path  is  greater  than  that  of  an  (A+a) -edged 
path,  where  a  is  any  positive  integer;  i.e.,  a  ^  1. 

Proof: 

Equation  (3.5)  can  be  written  for  an  (A+a) -edged  path. 
e-A  a 

_,.*  (A+a+i)-(e-A-a-i)       /e-A-a]  (A+a)_(e-A-a)jL 
K(i)pv     'pv        =     V  o   /  P     P 

i  -  0  \       / 

'"t")  p(A+a+i)?  <e-A-a-i) + rr)  p(A+a+2>  ?(e-A-a-2)  +-  (4.D 


+  /.-A-.  \    p(e-D  .1  +   /e-A-a]  pe 
le-A-a-1  I  \   e-A-a 


Comparing  equation  (4.1)  to  equation  (3.5),  we  see  that  theorem  (4.2) 
can  be  proved  if  it  can  be  shown  that,  using  definition  (4.2),  the  value 
of  equation  (3.5)  is  greater  than  the  value  of  equation  (4.1).   This  can 
be  done  by  showing  that  each  coefficient  of  equation  (4.1)  is  less  than 
that  of  equation  (3.5)  for  the  corresponding  power  of  p_. 
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Consider  the  general  term,  i  -  j.   If 


e-A-a  \  n(e-j)-j      ^  [    e-A    n(e-j)-j 

e-A-a-j    P     FJ      <       a-A-j   P     ^  (4.2) 


then 

e-A-a  \  /  e-A 

<  (         )  (43) 

e-A-a-j  /  \  e-A-j 


then 


(e-A-a).' s  <4-4) 

(e-A-a-j)  .'(e-A-a(e-A-a-j)).'  ^v 


(e-A).' 


(e-A^j)T(e^A(e^A-j))T 


then 


(e-A-a) (e-A-a-l)(e-A-a-2).  .  .  (e-A-a-j  1) 

j 


< 


(e-A)(e-A-l)(e-A-2)  .  .  .  (e-A-j  1) 

j! 


(4.5) 


There  will  be  j  number  of  terms  in  the  numberator  of  both  sides  of 
equation  (4.5),  each  term  of  which  is  greater  than  one. 
Since, 

e-A-a    ^    e-A 

e-A-a-1  <    e-A-1 

e-A-a-2  <    e-A-2 


e-A-a-j+1  <   e-A-j+1 
Then  equation  (4.5)  is  true.   From  it  equations  (4.4),  (4.3),  and 
(4.2)  follows,  and  the  theorem  is  proved. 
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Example  4.5: 

The  following  tables  are  the  numerical  values  of  the  2-terminal 
reliability  function  discussed  in  examples  (3.1  and  2)   Let  the  branch 
reliabilities  be  all  equal  and  have  the  reliability  0.9-   Columns  3,4, 
and  5  are  the  values  when  one  branch  is  removed  from  the  graph.   Column 
6  is  the  value  of  the  2-terminal  reliability  function  if  the  specified 
branch  has  the  reduced  reliability  as  shown. 


with  all 

without 

without 

without 

with 

branches 

branch 

branch 

branch 

P  =0.5 

present 
.981 

#1 
.81 

#2 
.90 

#3 
.90 

X 

Sc.y 

.905 

v.. 

.981 

.90 

.81 

.90 

.905 

q 
nx,z 

.981 

.90 

.90 

.81 

.905 

The  value 

of  q_ 

.972 

The  ideas  of  section  (3.3)  can  be  used  to  find  the  value  of  the 
2-terminal  reliability  function  in  the  cases  of  examples  (3.1  and  2). 
Equation  (3.8)  is  applicable  in  all  the  cases  where  the  reliabilities 
of  the  branches  are  all  equal;  equation  (3.7)  where  they  are  not  all 
equal. 

In  all  of  these  examples  there  exists  one  1-edged  path  and  one 
2-edged  path.   This  means,  using  equation  (3.8), 
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2- 

r  r  = 


2- 

P  P 


.9 
(.9)(.9)(.l) 


The  value  of  q    is 
x,y 


.9 

.081 
.981 


Observe  that  this  is  the  same  result  obtained  in  example  (4.5). 

Using  equation  (3.7)  the  case  of  example  (3.1)  where  the  reliability 
of  branch  one  has  the  reduced  reliability,  p.  ™  0.5. 


V3rl 


.5 


p2p3*i  =  <-9><-9><-5> 


The  value  of  q    is 
x,y 


.5 

.405 
.905 


Observe  that  this  is  the  same  value  found  in  example  (4.5).  It  is  signi- 
ficant to  observe  that  in  this  example,  if  branch  two  had  the  reduced  re- 
liability, the  effect  on  the  value  of  q    would  not  be  as  much.   Consider 

x»y 

example  (3.1),  let  branch  two  assume  the  reduced  branch  reliability,  p  = 


.5,  and  calculate  the  value  of  the  2-terminal  reliability  function,  q 
using  equation  (3.7). 


x,y 


Q 
Q 

The  value  of 


r  r  r 
2  3  1 


=  P2P3Pl 


.9 
(.5)(.9)(.l) 


V 


.9 

.045 

.945 


Observe  that  this  value  is  higher  than  that  one  found  in  the  case  where 
branch  one  had  the  reduced  reliability:   .945  vice  .905. 
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5.   Application  of  Ideas  without  the  Aid  of  a  Computer 

Some  of  the  ideas  presented  in  the  preceeding  sections  can  be  ap- 
plied, without  the  benefit  of  a  computer,  in  the  analysis  of  a  communica- 
tion net  to  give  rough,  but  fast,  measurements  of  the  properties  of  a 
net.   The  following  subsections  discuss  important  properties  of  the 
communication  net  which  must  be  examined  in  the  event  a  change  in  the 
construction  of  the  net  is  required.   For  example,  the  commander  of 
several  communication  systems  may  order  that  one  link  be  removed  from  one 
system  and  be  given  to  a  second  system.  The  effect  on  the  system  caused 
by  the  removal  of  a  link  can  be  studied  qualitatively  without  the  aid  of 
a  computer. 

5.1  Branch  loading. 

If  the  graph  of  a  communication  net  is  simple  or  if  sufficient  time 
is  available,  all  the  paths  of  the  net  can  be  identified,  and  the  load- 
ing of  each  branch  can  be  found.   Immediately  this  indicates  that  the 
branch  which  is  used  most  often  should  not  be  given  up,  because  all  the 
paths  using  that  branch  are  also  given  up.   On  the  other  hand,  the  light- 
est loaded  branch  may  not  be  sacrificed  in  every  case  without  additional 
examination  of  the  topology  of  the  graph.   That  branch  could  be  the  only 
connection  to  a  terminal,  and  the  destruction  of  this  branch  would  iso- 
late the  terminal.   However,  in  certain  conceivable  situations,  the  price 
of  isolation  of  a  terminal  may  be  acceptable,  for  example,  security 
reasons.   The  full  effect  of  branch  loading  on  the  decision  to  remove  a 
branch  must  be  tempered  by  additional  factors. 

5.2  Path  Length. 

Theorem  (4.2)  proves  that  long  paths  are  not  as  reliable  as  shorter 
ones.   Therefore,  one  must  place  a  figure  of  merit  on  the  communication 
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path  in  relation  to  its  length.  The  removal  of  either  one  of  two  branches 
may  still  enable  communication  between  stations;  however,  in  one  case 
the  remaining  paths  are  longer  than  in  the  other.   Then,  the  choice  would 
be  to  remove  the  one  that  would  leave  the  shorter  path.   This  choice  may 
be  a  poor  one  in  light  of  the  other  factors  which  must  be  considered. 

5.3  Branch  Reliability. 

Naturally  the  designer  of  a  communication  net  will  want  to  avoid  the 
use  of  a  branch  with  low  reliability.   And,  if  one  link  has  to  be  removed, 
then  that  one  of  lowest  reliability  should  be  selected.   However,  this  is 
not  always  true,  especially  when  that  branch  is  the  only  link  to  a  cer- 
tain terminal.  Again  the  price  of  isolation  of  that  terminal  may  be  ac- 
ceptable in  a  situation  where  security  considerations  are  paramount. 
Likewise,  the  link  of  highest  reliability  may  be  chosen  to  be  removed 
from  the  given  net  for  use  elsewhere,  especially  in  cases  where  the 
branch  loading  is  low  and  there  exist  other  paths  of  communication  of  ac- 
c eptable  reliability. 

5.4  Secondary  Path  Terms. 

The  importance  of  the  secondary  path  product  terms  has  been  stressed 
before.   It  is  pointed  out  in  this  section  that  since  they  constitute 
the  majority  of  terms  in  the  2-terminal  reliability  function,  and  they 
are  inversely  proportional  in  number  to  the  length  of  path  generating 
them,  the  loss  of  a  short  communication  path  is  very  costly  in  terms  of 
reliability. 

The  number  of  k-edged  paths  which  will  be  lost  when  a  given  branch, 
b  ,  is  removed  can  be  predicted  in  certain  cases. 
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Theorem  5.1: 


When  branch,  b  ,  is  removed  from  a  graph,  one  1-edged  path  will  be 


removed . 


Proof: 


There  can  be  only  one  1-edged  path  between  terminals  x  and  y_  con- 
sisting of  the  branch  b.  alone.  Removal  of  the  branch,  b  ,  results  in 
the  removal  of  that  path. 


Theorem  5.2: 

Let  D  and  D  be  the  degree  of  terminals  x  and  y_  respectively. 
x      y 

Let  b.  be  a  branch  connecting  terminals  x  and  y_.  Let  c  be  the  number  of 
additional  branches  connecting  terminals  x  and  y_.   Then  the  number  of  2- 
edged  paths  between  any  pair  of  terminals,  but  using  branch  b.,  that  will 
be  removed  if  branch  b.  is  removed,  is  given  by 

N0  =  (D  +  D  -  2c  -2)  (5.1) 

2     x    y 

Proof: 

The  degree  of  a  terminal  is  the  number  of  edges  incident  to  it. 

A  2-edged  path  using  the  branch,  b.,  must  either  start  at  terminal  x 

and  proceed  through  y_,  or  start  at  terminal  y_  and  proceed  through  x. 

In  the  former  case,  there  are  (D  -1)  number  of  possible  2-edged  paths 

minus  the  c   number  which  would  form  a  circuit,  since  they  are  connected 

to  terminal  x  also.   Similarly  in  the  latter  case,  there  are  [(D  -l)-c] 

number  possible.   The  sum  is  N_: 

N0=D  -  1  -  c  +  D  -1-c 
2    y  x 


N  -  D  +  D   -  2c  -2 
2    y    x 
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Example  5.1 


D  -  5 
x 

D  -  5 

y 

c  =  1,  namely  b_ 


Figure  5.1 


N2  =  5  +  5  -  2(1)  -2  =  6 


Namely, 

path  sets  Jb^      ,     [blb5]     ,    [bjb^     ,     Jb^J    ,     Jb^J     , 

s„d     Jbt       b8^ 

Rule  5.1. 

A  rule  which  establishes  a  series  of  operations  that  indicates  the 

number  of  3-edged  paths  between  any  pair  of  terminals  containing  the 

branch,  b  ,  that  will  be  removed  if  the  branch,  b  ,  is  removed.   It  is 

stated  here  without  proof;  however  an  example  is  given. 

Step  #1   Starting  at  node  x,  proceed  to  the  terminal  of  a  branch 

incident  to  node  x. 

Step  #2   Call  this  node  #1  having  degree  D1 .   Let  c  be  the  number 

of  branches  having  terminals  x  and  y_,  but  not  being  branches, 

b . .   Let  d.  be  the  number  of  branches  having  node  y_  and 

node  i   as  terminals.  Then 

D.  +  D  -  2  -  c  -2d.  -  N„  (5.2) 

i    y  i    3 

Step  #3   Proceed  to  all  other  nodes  of  the  branches  incident  to 

node  x>  except  the  node  y_,  of  course.   Repeat  the  calcu- 
lation of  equation  (5.2),  adding  the  result  to  the  pre- 
vious value  of  N». 
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Step  #4 


Proceed  to  the  terminals  of  the  branches  incident 

to  node  y_  and  whose  terminal  nodes  are  not  those 

previously  numbered  or  node  x*   Count  the  number  of 

branches  which  are  not  incident  to  x  or  ^  or  the  nodes 

previously  numbered  in  steps  (2)  or  (3).  Add  this 

number  of  the  previous  value  of  N_.   N_  is  the 

number  of  3-edged  paths  that  will  be  removed  if 

branch,  b.,  is  removed, 
l 


Example  5.2; 


Figure  5.2 
c  =  1  D 


N, 


Dx  -  3 
3+6-2-1-0  =  6 


d,  - 


D2  =  3 


•    d2  =  1 


3  +  6  -  2  -  1  -  2(1)  =  4 


N, 


N  -  6  +  4  =  10 
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D3  -  2,       d3  -  1 
2  +  6  -  2  -  1  -  2(1)  =■  3 
N        H_  -  10  +3  -  13 

\  -   3    >    \-1 
3+6  -2-1-  2(1)  =  4 

N3        N3  =  13  +  4  -  J£ 

Next  node  5  where  2  branches  qualify 


N  =  17  +  2  =  19 


ANS:   19 

5.5  Number  of  Paths 

Theorem  (3.3)  proves  that  the  primary  path  term  will  not  be  a 
duplicate  of  any  entry  in  the  Table  of  Combinations  regardless  of  the 
number  of  paths  that  exist  between  terminals  x  and  y_.  Therefore,  each 
path  will  contribute  some  value,  in  the  sense  of  definition  (4.2),  to 
the  2-terminal  reliability  function,  q   ,  and  its  removal  will  deminish 

x»y 

this  value. 

An  examination  of  the  paths  between  terminals  x  and  y_  alone  is  not 
sufficient  when  n-terminal  reliability  is  concerned.  What  may  appear  to 
be  an  expendable  path  between  terminals  x  and  y_,  may  be  a  vital  path 
between  two  other  terminals. 

While  a  large  number  of  paths  between  the  terminals  of  considera- 
tion indicates  that  the  loss  of  one  path  would  not  cause  isolation,  the 
length  of  the  paths  may  be  such  as  to  result  in  low  reliability.  All 
the  factors  mentioned  previously  continue  to  apply  even  when  there  are  a 
large  number  of  paths.   This  is  also  true  when  there  are  a  small  number 
of  paths. 
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5.6  Addition  of  a  branch. 

The  addition  of  one  branch  to  a  graph  is  equivalent  to  the  addi- 
tion of  one  communication  link  to  a  net.  This  can  be  viewed  as  a  de- 
sign problem  in  which  there  are  eight  cases  to  be  studied: 

What  should  the  second  branch  terminal  be,  if  the  given  objective 
is  to: 

the  first  branch  terminal  is 
required  to  be: 


maximize 

q 
x,y 

(1) 

terminal 

X 

(2) 

terminal 

a,  not  x  °r  y_ 

(3) 

optional 

maximize 

\ 

(4) 

terminal 

a  in  set  k 

(5) 

terminal 

a  not  in  set  k 

(6) 

optional 

maximize 

\ 

(7) 

(8) 

terminal 
optional 

X 

These  eight  cases  can  lead  to  many  more  if  the  reliability  of  the  branch 
to  be  added  and  the  reliabilities  of  the  other  branches  are  permitted  to 
change.  The  following  discussion  of  the  above  cases  assumes  that  all 
branch  reliabilities  are  the  same,  and  equal  p. 

Case  (1) : 

Corollary  (4.1)  specifies  terminal  y_  as  the  second  branch  terminal. 
Case  (2): 

Find  all  the  paths  between  terminals  x  and  a,  and  y_  and  a.  Choose 
terminal  x  or  y_  depending  on  which  has  the  fewer  paths  to  terminal  a. 
This  will  provide  the  most  paths,  with  the  fewest  branches,  between  x 
and  y_. 
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Case  (3) : 

Corollary  (4.1)  specifies  the  selection  of  terminals  x  and  y_. 
Case  (4) : 

Find  all  the  paths  between  all  the  pairs  of  terminals  of  set  k.   If 
any  one  pair  has  an  unusually  small  number,  examine  the  topology  of  the 
graph  in  light  of  sections  (5.1-4).  Base  the  selection  on  the  unusual 
topology.  Do  the  same  thing  in  the  event  one  pair  has  an  unusually  large 
number  of  paths  between  them.  Otherwise,  locate  the  terminal  that  has 
the  most  number  of  paths  to  terminal  a.  Select  this  as  the  second  ter- 
minal for  the  branch  to  be  added.  This  choice  will  add  one  1-edged  path, 
but  also  the  maximum  number  of  2-  and  3 -edged  paths  between  the  other 
terminals  in  the  set  k.  This  will  increase  the  value  of  q,  the  most. 
Case  (5): 

Find  all  the  paths  between  all  the  pairs  of  terminals  of  set  k  and 
terminal  a.   Examine  the  topology  of  the  graph  in  light  of  sections  (5.1- 
5).   Especially,  look  for  vertices  of  degree  one  or  two.  Unusual  con- 
figurations must  be  solved  in  this  manner.  Otherwise,  select  as  the 
second  terminal  that  one  in  set  k  that  has  the  largest  number  of  paths  to 
terminal  a.  See  Case  (4) . 
Case  (6) : 

Find  all  the  paths  between  all  the  terminal  pairs  in  set  k.   Choose 
as  the  terminals  of  the  branch  to  be  added  those  between  which  there  al- 
ready exist  the  largest  number  of  paths.   This  will  not  be  true  if  an 
examination  of  the  topology  reveals  an  unusual  nature  to  the  net.   See 
Case  (4). 
Case  (7); 

Find  all  the  paths  between  all  the  terminals  of  the  graph  and  ter- 
minal x.   Examine  the  topology  of  the  graph.   If  nothing  extradorinary 
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appears,  choose  as  the  second  terminal  for  the  branch  to  be  added  that 
terminal  which  already  has  the  most  paths  to  terminal  x. 
Case  (8); 

Find  all  the  paths  between  all  the  terminal  pairs  of  the  graph. 
Study  the  topology  of  the  net.   Choose  the  terminal  pair  between  which 
there  already  exist  the  largest  number  of  paths. 

If  there  exist  a  large  number  of  paths  between  two  terminals  of  a 
graph,  the  paths  are  usually  of  long  length.  This  means  that  they  rely 
on  many  other  branches,  and  so  are  of  low  reliability.   The  adding  of  a 
1-edged  path  between  these  terminals  immediately  increases  the  2-termin- 
al  reliability  between  them,  and  also  provides  an  increase  in  the  number 
of  2-edged  paths  between  other  pairs  of  terminals.   In  fact,  it  increases 
the  number  of  all  k-edged  paths.   Since  there  were  all  those  paths  in 
existence,  the  accessibility  of  the  new  branch  by  the  other  terminals 
produces  a  maximum  number  of  paths  between  all  terminals. 

Lemma  (3.1)  shows  that  an  isolated  terminals  I_,  will  cause  both 
the  2-  and  the  n-terminal  reliability  function  to  be  zero.   When  an 
isolated  terminal  exists,  clearly  the  topology  of  the  graph  is  unusual. 
The  isolated  terminal  must  be  selected  as  the  second  terminal  for  the 
branch  to  be  added  in  each  case  that  the  terminal  I  is  in  the  set  k. 
This  statement  covers  ail  the  cases  since  the  set  k  includes  all  situa- 
tions from  2  to  n  number  of  terminals. 

Example  5.3: 

Three  illustrations  of  a  2-terminal  reliability  function  are 

discussed  for  the  graph  of  figure  (5,3):  q  ,,  q,   ,  q   .   All  the 

31   b  >8   c  »8 

branch  reliabilities  are  equal,  p  ■  .9. 
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Figure  5.3 


qa,f 

qb,g 

qc,g 

Branch  Loading: 


with  all 
branches 
present 

.9859 

.9940 

.9963 


without  without    without 
branches  branches   branches 

#3      #6 


.9685 
.9621 
.9854 


.9755 
.9922 
.9937 


#9 
.9750 
.9897 
.9830 


qa,f 
qa,g 
qc,g 


Branch  number:      123456789  10 

Branch  load:          6444435344  10  paths 

1343333431  8  paths 

3434343333  9  paths 


Total  loading:  10  11  11  11  10  10  11  10  10  8 
In  the  illustrations  given,  branch  #3  is  not  the  one  to  be  removed. 
However,  branch  #10  should  not  be  removed  either  because  this  would  cause 
node  A  to  rely  completely  on  branch  #1  for  communication.   The  low  level 
of  branch  loading  for  branch  #10  indicated  that  it  is  ineffectively 
connected,  not  that  it  can  be  sacrificed. 
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6.   Digital  Computer  Program. 

6.1  Program  Specifications. 

A  digital  computer  program  has  been  written  in  the  FORTRAN  IV 

Language  suitable  for  use  on  the  IBM  System  360  computer  complex.   The 

primary  purpose  of  the  program  is  to  compute  the  2-terminal  reliability 

function,  q   .   This  program  is  included  as  Appendix  I  .   The  func- 
T<,y 

tional  block  diagram,  figures  (6.1)  (6.2)  and  (6.3)  is  fully  explained 
in  section  (6.2) . 

The  capabilities  of  the  program  go  beyond  the  computation  of  the  2- 
terminal  reliability  function,  q   .   After  finding  q   ,  the  program 

^»y  x»y 

can  calculate  its  numerical  value;  then,  it  can  recalculate  it  under 
the  condition  that  a  branch  of  the  graph  has  been  removed,  and  that  a 
branch  has  a  different  value  of  reliability.   (See  options,  Section 
6.1.1.)  After  the  second  and  any  succeeding  runs,  (a  run  is  defined  as 
the  computation  of  one  2-terminal  reliability  function, )q,,  where  k  is 
the  set  of  terminal  pairs  of  the  preceeding  runs,  can  be  computed,  and 
its  numerical  value  found.   In  finding  q   ,  all  path  sets  between 

x»y 

terminal  pairs  are  found;  this  information  is  given  as  a  preliminary  out- 
put. 

6.1.1  Input 

The  input  required  by  the  computer  program  consists  of  three  types: 

(1)  The  number  of  edges,  and  the  number  of  vertices  of  the  graph 
are  the  first  data  required.   The  number  of  fundamental  circuits  is  cal- 
culated from  this  data  using  the  relation: 

(e  -  n  +  1)  ■  w 
Also  on  the  first  data  card,  the  option  SEL1  is  provided  which,  if  used, 
enables  the  branch  reliabilities  to  be  input.   This  will  allow  the 
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calculation  of  the  numerical  value  of  q   ,  and  later,  the  numerical 
value  of  q,  ,  if  requested. 

(2)  The  w  number  of  fundamental  circuits  are  the  next  input 
data.   If  the  branch  reliabilities  are  to  be  used,  they  are  read  in  next, 

(3)  The  data  required  by  each  run  is  read  next.  This  data 
consists  of  three  main  parts: 

a)  the  terminal  vertices  (numbered) ; 

b)  one  path  between  these  vertices; 

c)  option  selections: 

1)  SEL  -  additional  runs  will  follow 

2)  SEL2  -  find  q 

3)  SEL4  -  calculate  the  numerical  value 

4)  SEL5  -  of  q    with  the  specified  branch 

5)  SEL6  -  having  been  removed  from  the  graph 

6)  SEL7  -  calculate  q    using  the  specified  value 

x»y 

of  branch  reliability  for  the  specified 

branch 
Precise  details  of  the  format  for  the  data  cards  are  provided  in  sec- 
tion (6.3). 

6.1.2  Output. 

The  output  provided  by  the  computer  program  is  fully  labelled  for 
easy  reading.  All  the  input  data  is  printed  with  the  generated  data  to 
provide  a  complete  record  of  the  run.   The  output  of  generated  data  in- 
clude: 

1)  Table  of  Combinations; 

2)  path  matrix; 

3)  branch  loading  information; 
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4)  number  of  k-edged  paths; 

5)  numerical  values  of  q    and  q,  . 

Mx,y     Mk' 

At  the  conclusion  of  the  computer  program,  a  summary  of  the  branch  load- 
ing and  the  number  of  k-edged  paths  is  printed.   A  message  explaining 
the  error  stop  is  printed  as  necessary. 

6.1.3  Theoretical  Limitations. 

In  its  present  form,  the  computer  program  is  limited  to  a  graph 
containing  no  more  than  15  edges,  or  10  nodes,  or  10  fundamental  cir- 
cuits.  There  must  be  at  least  one  circuit  in  the  graph.   Other  limita- 
tions include  a  maximum  of  4000  entries  in  the  Table  of  Combinations  for 
all  runs,  and  a  maximum  of  210  runs.   The  reason  for  these  limits  is  the 
storage  requirement,  especially  that  for  the  Table  of  Combinations.  Dif- 
ferent programming  techniques  might  be  able  to  increase  these  limits 
slightly.   The  fact  remains  that  the  powers  of  2  increase  very  rapidly, 
regardless  of  the  programming  technique. 

6.2  Block  Diagram. 

The  functional  block  diagram,  figures  (6.1,2  and  3),  is  presented 
to  give  a  better  description  of  the  operation  of  the  computer  program. 
The  blocks  are  identified  in  the  listing  of  the  computer  program,  Ap- 
pendix I. 
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/StartN 


Initialize 
Computer 


Convert  to 
Logical  bytes 


Print 
out 


Read  Branch 
reliabilities 


Initialize  con- 
stants for  run 


I 


Initialize  con- 
stants (rerun) 


(Read  input  for"\ 


Convert  to 
logical  bytes 


10 


I 


Build  combina- 
tion matrix, 
BFC 


11 


Ringsum 
PFC  ©   BFC 


6 


Figure  6.1 
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^9 


Write  original 
path  in  PFC 


13 


Generate  the 
path  product 
terms 


Ik ± 


Print  out 
path 
product 
terms 


Generate  path 
product  terms 


Printout 

path 

products 

terms 


Figure  6.2 
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19 


Printout 
path  set 
matrix 


Figure  6.3 
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Description  of  Block  Diagram 

Block  0:   Control  cards  and  dimension  and  format  statements  are  required 

to  prepare  the  computer  for  running. 
Block  1:   Input  data  are  read  including  number  of  branches  and  nodes, 

and  the  fundamental  circuits. 
Block  2:   Graph  size  is  tested  to  be  within  program  limits. 
Block  3:  Input  data  is  converted  to  logical  elements  to  facilitate 

faster  calculation. 
Block  4:   If  branch  reliabilities  are  to  be  used,  they  are  read  in 

now  into  DREL  matrix. 
Block  5:   Preliminary  printouts  are  made  for  more  complete  data  refer- 
ence for  the  final  output. 
Block  6:  Values  required  for  first  run  of  the  program  are  set. 
Block  7:  Values  required  for  succeeding  runs  are  set. 
Block  8;  Data  required  for  the  run  is  read  including  terminals  and 

path. 
Block  9:  Path  data  is  converted  to  logical  elements,  and  a  storage 

area  prepared  for  later  use. 
Block  10:  All  possible  combinations  of  the  given  circuits  are  generated 

and  stored  in  matrix  BFC. 
Block  11;  The  ringsum  operation  is  performed  on  the  given  path  and 

each  entry  in  matrix  BFC,  and  the  result  is  written  in  BRP. 
Block  12;  The  original  path  is  written  into  the  path  matrix,  PFC. 
Block  13;  The  primary  and  all  the  secondary  path  product  terms  are 

generated  into  OPC. 
Block  14;  The  path  product  terms  are  printed  out. 
Block  15:  Matrix  BRP  is  searched  for  pure  path  sets. 
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Block  16:  The  primary  and  secondary  path  product  terms  are  generated 

for  the  newly  formed  path  set,  and  written  into  QTEMP. 
Block  17:  Each  entry  in  QTEMP  is  compared  to  all  entries  in  QPC. 

If  a  duplicate  is  not  found,  the  set  from  QTEMP  is  written 

in  QPC. 
Block  18:  Those  sets  added  to  QPC  are  printed  out.  The  number  of 

entries  in  QPC  and  the  number  in  PFC  are  printed. 
Block  19:  The  path  matrix,  PFC,  is  printed  out. 

Block  20:  The  various  OPTIONS  are  examined.   If  Q.  ,  or  Q  are  to  be 

i>  3  n 

found, transfer  is  made;  if  another  run  it  to  be  made,  trans- 
fer is  accomplished. 

Block  21:  Before  stopping,  a  summary  is  printed  out  including  branch 

loading,  number  of  paths  between  each  node,  and  path  length. 

Block  22:  Error  stops  are  provided. 

Block  23:  Generate  matrix  Q  from  path  product  terms  common  to  the  first 

n 

2  q,   matrices. 
i,v 

Block  24:  Generate  matrix  Q  from  the  previous  Q  and  the  newly  formed 

m  n 

q   .  for  this  run. 

Block  25:  Print  out  matrix  Q  . 

n 

Block  26:  If  the  branch  reliabilities  are  provided,  determine  the  value 

of  Q  .   Transfer  control  to  block  20  to  check  options, 
n  r 

Block  27:  If  the  branch  reliabilities  are  provided,  determine  the  value 

of  q.  ..   Examine  the  options  SEL4,  SEL5,  SEL6,  and  determine 

the  value  of  q.  .  without  the  specified  branch  J.   Examine 
i,j 

option  SEL7,  and  compute  the  value  of  q.  .  with  the  revised 
value  of  branch  reliability.   Transfer  control  to  block  20  to 
check  other  options. 
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6.3  Data  Card  Configuration. 

The  following  is  a  description  of  the  data  cards  required  for  the 
running  of  the  computer  program  as  described  in  section  (6) . 


Data  Card 

?ype  Number: 
1 


Columns : 


4  -  6 


10  -  72 


Description 

and  specifications: 

Number  of  edges,  right  justified; 

max:  15 

Number  of  nodes,  right  justified; 

max :  10 

If  0  or  Blank:   no  branch  reliabilities 

are  to  be  used;  if  nonzero:  branch 

reliabilities  are  to  be  read  and  used. 

Disregarded. 


15 


15 

-  72 

1 

-  5 

6 

-  13 

14 

-  21 

22 

-  29 

30 

-  37 

38 

-  45 

46 

-  72 

Fundamental  circuit  set: 
Up  to  15  digits,  0  or  1,   left 
justified,  representing  the  edges  in 
the  circuit; 

1  if  edge  present 

0  if  edge  not  present 

Disregarded. 
Blank 

Decimal  point  and  7  digits  to  represent 
the  branch  reliability  of  the  succes- 
sive branches  1,  2,  .  .  .  ,  e 
One  field  for  each  edge.  Up  to  five 
edges  per  card. 
Disregarded. 
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Description 
Type  Number;       Columns:       and  Specifications: 

4  1-3         Identifying  number  of  starting  node 

4-6         Identifying  number  of  terminal  node 

7-9         if  0,  no  more  cards  will  be  read; 

if  not  0,  additional  cards  read. 

10  -  12        if  nonzero,  q,  will  be  computed  using 

the  previously  found  q.  .,   from  all 

l ,  j  s 

previous  runs  as  the  set  k;  if  0, 
then  this  is  not  done. 
13  -  27         left  justified,  up  to  15  digits, 

0  to  1,  representing  the  edges  in  a 
path  between  the  subject  terminals; 

0  if  edge  not  present 

1  if  edge  present 

only  e  number  of  digits  processed 
28  -  30        an  integer,  less  than  or  equal  to  e, 
representing  an  edge  that  has  been 
removed  from  the  graph.   A  value  for 

q   .  is  found  assuming  that  edge  is 

* » J 

removed  from  the  graph. 

31  -  33        another  edge  number  processed  as  that 
on  in  the  preceeding  field. 

34  -  37        Same  as  above 

38  -  40        An  edge  number  whose  value  of  branch 
reliability  it  is  desired  to  change 
temporarily,  and  recompute  q   .  with 
the  new  value  of  reliability. 
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Type  Number: 


Columns: 


41  -  48 


Description 

and  Specifications: 

A  decimal  point  and  7  digits  represent- 

the  branch  reliability  of  the  edge 

specified  in  the  preceeding  field. 
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7.   Conclusions. 

This  work  shows  that  the  reliability  function  as  defined  is  an  ef- 
fective instrument  in  the  analysis  and  design  of  communication  nets.   The 
theorems  that  are  developed  can  be  applied  in  many  cases  without  the  use 
of  a  digital  computer  to  provide  direction  in  the  modification  of  a  net. 
The  computer  can  provide  fast  data  to  quantify  the  effect  of  any  change 
to  a  communication  net  with  respect  to  k-terminal  reliability,  as  the 
results  of  the  computer  study  indicate. 

The  computer  program  can  be  modified  by  specialization  to  facilitate 
sensitivity  studies  for  branch  removal  or  branch  addition.   In  so  doing, 
the  limitations  on  graph  construction  can  be  modified.  Moreover,  equa- 
tion (3.7)  can  be  used  in  a  computer  program  to  study  the  2-terminal 
reliability  function  alone. 
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APPENDIX  I 
The  digital  computer  program,  as  described  in  section  (6)  is 
presented  as  Appendix  I.   The  programming  language  is  FORTRAN  IV  and 
the  program  is  suitable  for  use  on  the  IBM  System  360  computer  complex. 
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0001 
0002 


C 

c 
c 


C003 
0304 


0005 
00G6 
0007 

0008 

0009 
0010 
0011 
0012 
0013 
0014 
0015 

0  016 

0017 
0018 
0019 

0020 
0021 
0022 
0023 
0024 
0025 
002b 
0027 

C}28 


C 
C 


BLOCK   NUMBER   0 

IMPL  IC  IT  L0G1(.AL*HA-C,0-Q)  , 
1INTEGEP*4(E-N,P-Z) 

0  1  MENS  I  ON  X8F( 10,15) ,BC«(  10,15)  , 
13FC( lC24,15),BPP(1024,l5)f  PFCU024,  15)  , 
IXPFCI (210,15>TXY(210,15),ENO(210,15)t 
1GN0(210) ,QN0(4C0,15) ,QN1 (400,15) , 
1QTS( 15) ,0PC(40C0,15) , QTEMP ( 400 , 1 5 ) , 
1UN0121C)  ,  URUN1213)  ,  XPTHU5.)  , 
1PTSI  15)  , ST (2 10) , SP(210) , OR  EL  I  15) 


45 
46 


76 
77 

73 

79 

80 
^1 
82 
83 
84 
8  5 
86 

37 

8  8 
39 
90 
91 
92 
93 
94 
95 
96 
97 
98 

99 


FORMAT 

115H  WI 
FORMAT 

114H  WI 

125H  HA 
FORMAT 
FORMAT 
FORMAT 

133H  K- 
FflRMAT 

120HHS 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

UOH  LI 
FORMAT 

125H7 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

16H  FP 
FORMAT 


</,26H 
THOUT  BR 
(/,2  6H 
TH  BRANC 
VING  RFL 
(/,6H   S 
(/,22H 
(/,25H 
EDGED  PA 
< IH1,25H 
FOUND  IN 
(313) 
(1511) 
(216,151 
(5X,  1513 
(  1H1,2CH 
(5X,/) 
(  1H1.1CH 
NKS  AND1 
( /,25H 
8   9  1^ 
(/,29H 

(/,2CH 
(/,5H 
(/,20H 
(413,151 
(5X, 5F8. 
(/,30H 
(/,13H 
(/,11H 
(/,19H 
(  1H1, 8H 
0M1I 5,4H 
( /,10H  T 


THE  VALUE  OF  Q-SU8-( I ,J) 
ANCH1I3) 

THE  VALUE  OF  Q-SU8-( I ,J) 
H   113,/, 

IABILITY  VALUE1F8.6) 
TO  P 1 1 5 ) 

TOTAL  BRANCH  LOADING) 
IN  PAIRS,  THE  NUMBER  OF 
THS  AMD  BRANCH  LOADING) 
THE  TOTAL  NUMBFR  OF  PAT 
EACH  RUN) 


1,112) 
) 

THE  GIVEN  PATH  WAS) 

THERE  ARE  115, 
15, 6H  NUDES/) 
FDGE  12   3   4   5   6 
11  12  13  14  15/) 
THE  GIVEN  CIRCUIT  MATRIX 

THE  PATH  MATRIX  IS) 
P=1I3,5H    U=1I3) 
MATRIX  QPC  FOLLOWS) 
1 ,4I3,1F8.6) 
6) 

THF  BRANCH  RELIABILITIES 
Q-SUB-I,J  =1F10.7) 
C-SUB-K  = IF 10. 7) 
MATRIX  QK  FOLLOWS) 
PUN  NO.  116, 

T01  IS,/) 
HERE    ARE1I5,6H    PATHS) 


IS) 


ARE) 
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GO 

TO 

3 

GO 

TO 

3 

GO 

TO 

3 

GO 

TO 

■x 

GO 

TO 

3 

GO 

TO 

3 

c 

C  BLOCK   NUMBER   1,2,3,4,  AND  5 
C 

0029  RFAD(5,80)  E,N,SEL1 

0030  WPITE(6f86)  E,N 

0031  W=E-N+1 

0032  IF(E.GT.15) 

0033  IF(E.Lfc,  0) 

0034  IF(N.GT.IO) 

0035  IF(N.LE.  0) 

0036  IF(W.GT.IO) 

0037  IF(W.LE.  0> 

0038  00  71  1=1, W 

0039  71  READ(5,81) ( X8F ( I , J) ,J=1,E) 

0040  00  75  1  =  1, W 

0041  00  70  J=lrE 

0042  IF  (XBF( I , J) .FC.O)   BCM ( I , J ) =. F ALSE 

0043  IF(XBF( I ,J).E0.1 )   BCM( I , J ) = . TRUE . 

0044  70  CONTINUE 

0045  75  CONTINUE 

0046  WRITE(6,88) 

0047  00  72  1=1, W 

0048  72  WRITE(6,83) (XBF (I , J) , J=1,F) 

0049  WRITE  (6,85) 

0050  IF(SELl.EO.O)  GO  TO  10 

0051  RFAD(5,93MDRFL<  J)  ,J=1,F) 

0052  WPITE(6,94) 

0053  WRITE(6,93)(DPFL( J), J=1,E  ) 
C 

C  BLOCK   NUMBER   6  AND  7 
C 

0054  10  TIMS=0 

0055  OPT  =  .FALSE. 

0056  OTAG  =. FALSE. 

0057  QTAG1=. FALSE. 

0058  QTAG2=. FALSE. 

0059  GR=0 

0060  R=0 

0061  U=0 

0062  1  R=R+1 

0063  URAN=U+1 

0064  U=0 

0065  UV=C 
C 
C  BLOCK   NUMBER   8  AND  9 


C 
0066 


READ(5,92>ST(R),SP(R),SEL,SEL2,XPTH, 
1    SEL4,SEL5,SEL6,SEL7,0CHG 
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0067  WRITE(6,84) 

0068  WRITE(6t87) 

0069  WRITE(6,83J  XPTH 

0070  WRITE  (6,85) 

0071  00  54  1=1, E 

0072  IF  (XPTH( II.EQ.O)   PTS( I ) =. FALSE. 

0073  IF  (XPTH(I).EQ.l)   PTS ( I )= .TRUE. 
00  74  XY(RrtM=0 

0075  54  ENO(R,J)=0 
C 

C  3 LOCK   NUMBER   10 
C 

0076  IF(QPT)GQ  TO  51 

0077  OPT  =  .TRUE. 

0078  00  55  1=1, E 

0079  BFCUt  n*.FALSF. 

0080  55  BFC(2t I)-BCM(1VI) 

0081  K=2 

0082  G=3 

0083  IF(W.EQ.l)    GO  TO  51 

0084  DO  58  1=2, W 

0085  00  57  J=1,K 

0086  DO  56  F=1,E 

0087  A=BFC(J,F) 
0C88  3=BCM(I,F) 
0089  C  =. FALSE. 

C090  IF  (A.OR.BI   C=.TRUF. 

0091  IF  (A.AND.B)  C=. FALSE. 

0092  BFC(G,F)  =  C 

0093  56  CONTINUE 
C094  G=G«-1 

0095  57  CONTINUE 

0096  58  K=G-1 
C 

C  3LOCK   NUMBER   11 
C 

0097  51  T=GR+1 

0098  DO  38  1=1, K 

0099  DO  37  J=1,E 

0100  A=BFCU,J> 

0101  52  B=PTS(J) 

0102  C  =. FALSE. 

0103  IF  (A.OR.B)   C=.TRUE. 

0104  IF  (A.AND.B)  C=. FALSE. 

0105  37  3RP(I,J)=C 

0106  38  CONTINUE 
C 

C  BLOCK   NUMBER   12,  13,  AND  14 
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0107 
0108 
0109 
0110 
0111 
0112 
0113 
0114 
0115 
Ollt, 
0117 
0118 
0119 
0120 
0121 
0122 
0123 
0124 
0125 
0126 
0127 
0128 
0129 
0130 
C131 
0132 
0133 
0134 
0135 
0136 
0137 
0138 
0139 
0140 
0141 


014? 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 


c 

c 

c 


U=URAN 

UV=U 

DO  10  3  J=1,E 

QTS( J)=PTS(J) 

gPC(U»J)=PTS( J) 
103  PFC(T,J)=PTS( J) 

i)0  115  J=1,E 

IF(QTSCJ).)  GO  TO  115 

JTS( J)=.TRUE. 

USTOP=UV 

00  102  V=U»USTCP  , 

UV=UV+1 

00  101  JJ=1,E 

C=. FALSE. 

A=QPC< V,JJ) 

3=QTS(JJ) 

[F(A.OR.B)  C=.T*UE. 

OPC(UVTJJ)=C 
1C1  CONTINUE 
10  2  CONTINUE 

:JTS(  J)=.FALSF. 
115  CONTINUE 

U=UV 

WPITE(6,91) 

WRITE<6,87) 

H  =  0 

00  225  HO=UPANtUV 

H  =  H+1 

DO  224  J=l ,E 

XPFCKH,  J)  =  l 

IF(QPC<HO,J) )  GO  TO  224 

XPFCI (H,J)=0 

224  CONTINUE 

WRITE (6, 83) <XPFCI(H,J),J=1,E) 

225  CONTINUE 

3L0CK   NUMBER   15 

DO  228  I=1,K 

AD=. FALSE. 

DO  226  J=1,E 

IF(PTS( J) )  GO  TO  230 

226  CONTINUE 

IF(AD)  GO  TO  228 
GO  TO  236 
230  IF(8RP(  I,J)  )  GO  TO  234 
GO  TO  236 
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0151  234  AD=.TRUE. 

0152  GO  TO  226 

C   CIRCUIT  EXISTENCE  TEST 

0153  236  AD=. FALSE. 

0154  00  299  L=1,K 

0155  AD=. FALSE. 

0156  00  295  J=1,E 

0157  IF(8FC<LtJU  GO  TO  298 

0158  295  CONTINUE 

0159  IF  (AD)   GO  TO  223 

0160  GO  TO  300 

0161  238  AD=.TRUE. 

0162  GO  TO  295 

0163  298  IF  (BRP(ItJ)l  GO  TO  238 

0164  '  300  CONTINUE 

0165  299  CONTINUE 
C 

C  BLOCK   NUMBER   16,  17  AND  18 
C 

0166  T=T+1 

0167  V=l 

0168  DO  104  J  =  1,E 

0169  QTS( J)=BRP( I  ,J) 

0170  3TEMP< 1, J)=BRP<  I , J  ) 

0171  104  PFC(T,J)=BRP(  I  ,J) 

0172  00  108  J=1,E 

0173  IF(QTS( J) )  GO  TO  108 

0174  QTS( J)=.TRUE. 

0175  STEP  =V 

0176  DO  106  IV=1,STFP 

0177  V=V+1 

0178  DO  105  JJ=1,E 

0179  C=. FALSE. 

0180  A=QTEMP( IV, JJ) 

0181  B=QTS(JJ) 

0182  IF(A.OR.B)   C=.TRUE. 

0183  QTEMP(V,JJ)=C 

0184  105  CONTINUE 

0185  106  CONTINUE 

0186  QTSU>  =  .FAL5E. 

0187  108  CONTINUE 

0188  KTR=0 

0189  DO  113  H=1,V 

0190  DO  111  IV=URAN,U 

0191  DO  110  J=1,E 

0192  A=QTEMP(H,J) 

0193  d=QPC(IV,J) 

0194  IF(A.AMO.B)  GO  TO  110 
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0195 

0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0208 
02C9 
0210 
0211 
0212 
0213 
0214 
0215 


0216 
0217 
0218 
0219 
G22C 
0221 
0222 
0223 
0224 
0225 
0226 
C227 
0228 
0229 
0230 
0  231 
0232 
0233 
0234 
0235 
0236 
0237 


110 
111 


112 
113 


228 


126 


127 


IF(A.OR.B)    GCi  TO  111 

CONTINUE 

30  TO  113 

CONTINUE 

KTR=KTR+1 

UV=UV-H 

00  112  J=1,E 

XPFCI lltJ)=l 

IF(QTEMP(H,J) )  GO  TO  112 

XPFCI <  1,J)=0 

3PCUJV,  J)  =  QTEMP(H,JI 

WRITE (6,83)  (XPFCK  It  J)  ,J=1,E) 

CONTINUE 

u=uv 

IF(U.GT.4000)  GO  TO  5 

CONTINUE 

TNO=T-GR 

WPITE(6,85) 

WRITE (6,99)  TNG 

WRITE(6,89) 

WPITE(6,87) 

BLOCK   NUMBER   19 


TSTR  = 
H  =  0 

00  127 
^  =  H+1 

x=o 

00  126 


GR  +  1 

I  =TSTR,T 


J=l 


TO  128 


,E 
XPFCKHtJ)  =1 
IF(PFC<  ItJ)  )  GO 
XPFCI (H,J>  =  0 
CONTINUE 

WRITE  <6,83>  (  XPFCKH,  J),  J=l,E> 
XY(R?X)  =XY(R,X)+1 
CONTINUE 
WRITF(6,78) 
WRITE(6,87) 
WRITE(6,83> (XY(R,X) 
WRITE(6,83) (ENO(R, J) , 
WRITE  (6,85) 
URUN(R)=U 
UN0(R)=U-URAN+1 
GNO(R)=TNO 
GR  =  T 

BLOCK   NUMBER   20,  21  ANO  22 


X=1,E) 
J=1,E> 
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023Q 
0239 
0240 
0241 
0242 
0243 
0244 
0245 
0246 
0247 
0248 
0249 
0250 
0251 
0252 
0253 
0254 
0255 
0256 
0257 
0258 
0259 
0260 
0261 
0262 
0263 
0264 
0265 
0266 
0267 


0268 
0269 
0270 
0271 
0272 
0273 
0274 
0275 
C276 
0277 
0278 
0279 
028C 
0281 


C 

c 

c 


119 
121 
128 


53 


120 


IF 
IF 

0T 
IF 

GU 
EN 
X  = 
GO 
WR 
WR 
W 
DO 
UN 
00 
WR 

no 

UN 
WR 
WR 
WR 
WR 
WR 
ST 
RE 
GO 
RE 
GO 
RE 
WR 
GO 


(SEL1 

COTAG 
AG=.T 
(SEL 
TO  1 
0(R,  J 
X+l 

TO  1 
ITE(6 
ITEI6 
RITE( 

6  J  = 
0(  J)  = 

2  1  = 
ITE<6 

7  J  = 
0U)  = 
ITE(6 
ITE(6 
ITE(6 
ITE(6 
ITE(6 
OP 
STP  =  9 

TO  8 
STP  =  9 

TO  8 
S1P  =  8 
ITE(6 

TO  9 


.NE.O)  GO  TO  15 

)   GO  TO  120 

RUE. 

.EQ.O)   GO  TO  53 

)=ENG(R,J)+1 

26 

,79) 

,83MGNC(  H  ,  I  =  1,RJ 

6,78) 

1,E 

0 

l.R 

,98)  I,ST< I  ),SP< I) 

1,E 

UNO! J)+ENO( I ,J) 

,83)(XY( I, J),J=1,E) 

,83) (ENO( I, J)  , J=1,E) 

,85) 

,77) 

,83) <UNC( J)  ,J  =  1,E) 

9 

99 

Q  U 

,76)  PESTP 


BLOCK   NUMBER   23 


IF  (SEL?. EQ.O)   GO  TO 

IF(OTAGl)   GO  TO  150 

QTAG1=.TRUE. 

Z  =  l 

7.A1  =  UN0(1)   • 

ZA?«ZAi+l 

ZA3=ZA1+UN0(2) 

00  125   YZ=1,ZA1 

00  124  YY=ZA2,ZA3 

00  123   J=1,F 

A=QPC(YZ,J) 

4=QPC(YY,J) 

IF  (A.  AND.  B)   GO  TO  12?. 

IF  (A.OR.B)   GC  TO  124 


121 
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0282 

1?2 

QN0( Z,J)=A 

0283 

123 

CONTINUE 

0284 

7  =  14-1 

0285 

124 

CONTINUE 

0286 

125 

CONTINUE 

0287 

YR  =  2 

0288 

C 

C 

c 

150 

GO  TO  161 

3L0CK   NUMBER   ?4  AND  25 

0289 

7STR=URUN( YR)+1 

0290 

YR=YR+1 

0291 

ZMX=Z-1 

0292 

ZAR=URUN(YR) 

0293 

Z  =  l 

0294 

00  160  YZ=ZSTRTZAR 

0295 

00  159  YY=1,ZMX 

029o 

DO  158   J=1,E 

0297 

A=QPC<YZ,J> 

0298 

IF(QTAG2>   GO  TO  154 

0299 

3=QNC(YY,J) 

0300 

151 

IF(A.AND.B)   GO  TO  152 

0301 

IF(A.OR.B)   GO  TO  159 

0302 

152 

IF1QTAG2)   GO  TO  156 

0303 

QN1(Z»J)=A 

0304 

GO  TO  158 

0305 

154 

3=QN1( YY,J) 

0306 

GO  TO  151 

0307 

156 

QNO< Z,J)=A 

0308 

GG  TO  158 

0309 

158 

CONTINUE 

0310 

Z  =  Z+1 

0311 

159 

CONTINUE 

0312 

160 

CONTINUE 

0313 

QTAG2=.NOT.QTAG2 

0314 

161 

IF(YR.EQ.R)  GO  TO  162 

0315 

GO  TO  150 

0316 

16? 

ZMX=Z-1 

0317 

WPITE(6,97) 

0318 

DO  27  YD=ltZMX 

0319 

DO  26  J=i,E 
XPFCI (It J)=l 

0320 

0321 

IF(QTAG2>  GO  TO  24 

0322 

TF(QN0(YD,J) )  CO  TO  26 

0323 

25 

XPFCI ( 1,J)=0 

0324 

26 

CONTINUE 

0325 

WRITE  (6,  83 )  (XPFCM  1,J)  ,  J*l 

tE) 

0326 

27 

CONTINUE 
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0327  GO  TO  164 

0328  24  IF(QNKYDtJ))  GO  TO  26 

0329  GO  TO  25 
C 

C  BLOCK   NUMBER   26 
C 

0330  164  IF(SEL2.NE.2)  GO  TO  121 

0331  ZMX=Z-1 

0332  DSUM=0.0 

0333  HO  168  Y0=1,ZMX 

0334  0TEMP=1.0 

0335  00  167  J=1,E 

033d  IF(0TAG2)  GO  TD  165 

0337  IF(QNO(YDf J> )     GOTO  163 

0338  166  UEMP=DTEMP*(1.0-DREL< J) ) 
0  3  39-  167  CONTINUE 

0340  DSUM=DSUM+DTEMP 

0341  168  CONTINUE 

0342  WRITE(6,96)  DSUM 

0343  GO  TO  121 

0344  163  OTEMP=DTEMP*DREL(J) 
0  345  GO  TO  167 

0346  165  IF(QNl(YOrJ) )  GO  TO  163 

0347  GO  TO  166 
C 

C  BLOCK   NUMBER   27 
C 

0348  14  DTEMP=DTEMP*OREL(J) 

0349  GO  TO  16 

0350  15  DSUM=0.0 

0351  USTP=URUN(R) 

0352  00  17  UO=URANTUSTP 

0353  0TEMP=1.0 

0354  00  16  J=1,E 

0355  IF(QPC(UD,J) )  CO  TO  14 

0356  OTEMP=DTEMP*( 1.0-DREU J)  ) 

0357  16  CONTINUE 

0358  OSUM=DSUM+DTEMP 

0359  17  CONTINUE 

0360  WPITE<6,98>  R,ST(P),  SP<R> 

0361  WRITE<6,95)  DSUM 

0362  WRITE(6,85) 

0363  IF(SFL4.LE.0)   GO  TO  31 

0364  SELX=SEL4 

0365  30  IF(SELX.GT.15)   GO  TO  3 

0366  0TEM1=DREL(SELX) 

0367  DREL(SELX»=C.O 

0368  OSUM=0.0 


73 


0369  30    32    UD=URAN,USTP 

0370  IF(QPC(U0,SELX)  )    GO    TQ    32 
0  3  71                                       OTEffP=1.0 

0372  30    33    J=1,E 

0373  IF(QPC(UD,J) )     CO    TO    36 

0  3  74  nFMP=DTFMP*(  1.0-DR  EU  J)  ) 

0375  33    CONTINUE 

0376  0SUM=DSUM*OTFMf 

0377  32    CONTINUE 

0378  a'RITE(6,45)     SELX 
037^                                        'WRITE  <6t95)    OSUM 

0380  wRITE(6t85) 

0381  IRELi SFLX)=DTF^1 

0  382  31     IF(SEL5.LE.O)       GO    TO    34 

0383  SfLX  =  SF|_5 

03ti4  Sf-L5=0 

03  85  GO  TO  30 

0386  36     1Tf-MP  =  HTFMP*DRFL(  J  ) 

0  387  GO    TO    3  3 

0388  34     IF(SFL6.LE.O)       GO    TO     35 

0389  StLX=SEL6 
039C  SEL6=0 
0391                                         GO    TO    30 

C392       —  35     IF<SFL7.LE.O)       GO    TO    119 

3393  IF(SFL7.GT.15)        GO    TO    3 

0394  0TEM1=DREL(SEL7) 

0395  OREL* SHL7)=DCHG 

0396  0SUM=0.0 

0397  20  44  UD=URAN,USTP 

0398  3TFMP=1.0 

0  3  99  00  43  J=1,F 

04-00  IF<OPCUiD»J))  GO  TO  42 

0401  )TFyp=[)TFMP*(  1  .^-JREL( J)  ) 

0402  ^3  CONTINUE 

0403  -)SUM  =  OSUMfDTFMP 
04^4  44  CONTINUE 

0405  <*«ITF(6,46)  SE17,DCHG 

04C6  .'RITE  (6,  95)  DSUM 

04^7  WPITF(*,85) 

0408  >)R6USEL7)-DTEM1 

0409  '30    TO     119 

0410  <*1    DTEMP=DTEMP*DRFU  J  ) 

0411  n,D    TO    43 
C 

C 
^412  END 
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C  GR  =  TOTAL  NO  OF  PATHS  AT  START  OF  RUN 

C  T  IS  USEO  AS  A  RUNNING  INDEX 

C  GNO(R)  IS  USEO  AS  A  TOTAL  FOR  THE  RUN 

C  IF  SEL  IS  ZEROt  COMPUTER  WILL  HALT  AFTER 

C  THIS  RUN;   IF  NONZERO,  THEN  AN  ADDITIONAL 

C  DATA  CARD  WILL  BE  READ  AFTER  THIS  RUN 

C  IF  SFL2  IS  ZERO,  NO  Q-SUB-K  FUNCTION 

C  WILL  BE  FOUND;   IF  NONZERO,  THEN  Q-SU6-N 

C  WILL  BE  DETERMINED 

C  NOW  COMPUTE  Q-SUB-N 

C  UNO(R)  CONTAINS  LENGHT  OF  Q-SUB-I,J(P)   ' 

C  R  IS  THE  RUN  NUMBER 

C  ST  IS  THE  STARTING  TERMINAL 

C  SP  IS  THE  ENDING  TERMINAL 

C  W  IS  THE  NO.  OF  INDEPENDENT  CIRCUITS 

C  K  IS  THE  NO.  OF  POSSIBLE  COMBINATIONS 

C       OF  THF  GIVEN  CIRCUITS 

C  GR  IS  USED  AS  A  RUN-TOTAL-ER 

C  STO°  99  MEANS  THE  NUMBER  OF  NODFS,  EDGES,  OR  CIRCUITS 

C  IS  GREATtR  THAN  THE  LIMITS  OF  15,  10,  AND  10 

C  STQD  999  MEANS  THE  LIMIT  IF  Q-SUB-K  HAS  BEEN  REACHED 

C  STOP  888  MEANS  TOG  MANY  TERMS  HAVE  BEEN  PRODUCED 

C  AND  ALL  RESULTS  APE  NOW  IN  DOURT 

C  SEL  IS  THE  SELECTION  CODE 

C  IF  SEL1  =  0,  NO  bFANCH  RELIABILITIES  ARE  READ 

C  IF  SELl  NOT  ZERO,  THE  VALUE  FOR  Q-SUBI,J  IS  FOUND 

C  IF  SEL2  =  2,   THE  VALUE  FOR  Q-SU8-N  IS  FOUND 

C  MATRIX  XrSF  CONTAINS  THE  INDEPENDENT  CIRCUITS 

C  MATRIX  BCM  IS  XBF  AS  LOGICAL  ELEMENTS 

C  MATRIX  XPTH  CONTAINS  THE  GIVEN  PATH  SET 

C  MATRIX  BFC  CONTAINS  ALL  POSSIBLE  COMBIN- 

C  AT  IONS  OF  BCM 

C  MATRIX  BPP  CONTAINS  THE  RINGSUM  OF  THE 

C  PATH  AND  i^FC 

C  MATRIX  PFC  IS  THE  PATH  MATRIX,  INDEXED 

C  BY  T  OVER  ALL  RUNS 

C  MATRIX  XPFCI  IS  THE  INTEGER.  VERSIUN  OF  PFC 

C  MATRIX  XV  COUNTS  THE  NUMBER  OF  K-EOGED  PATHS 

C  MATRIX  ENO  CONTAINS  BRANCH  LOADING 

C  MATRIX  GNO  HAS  THE  NUMBER  OF  PATHS 

C  MATRIX 

C  THIS  ADDS  UNIQUE  PATH  PRODUCT  TERMS 

C  TO  THE  COLLECTION  NOW  IN  QPC 

C  BUILD  THE  PATH  PRODUCT  TERMS 

C  STORE  FCK  LATER  TESTING  OF  UNIQUENESS 
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APPENDIX  II 
The  problems  solved  in  examples  (3.1-3)  and  (4.5)  are  combined 
and  solved  by  the  digital  computer  program.   The  output  of  the  running 
of  the  program  is  included  as  Appendix  II. 
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THERE  ARE     3  LINKS  AND     i  NODFS 


THE  GIVEN  CIRCUIT  MATRIX  IS 
I   1   1 


THE  BRANCH  RELIABILITIES  ARE 
9  C  CO  C  0 Q  .  90  QOOnXL.  9O0C  C  0 


77 


THE  GIVEN  PATH  WAS 
E0GF~1   2   ?   4   5   6   7   8   9  10  11  12  13  14  15 

1   P   C   3   (TO   0   0   0   0   0   C   0 D   0 


MATRIX  QPC  FOLLOWS 


EDGE  I   2   3   4   5   6   7   8   9  1C  1 1  1 2  13  14  15 


1 

0 

■e 

1 

1 

o 

1 

0 

l 

1 

1 

l 

0 

1 

i 

THe*f"  ARE     2  PATHS 
THE  PATH  MATRIX  IS 


EDGE  I   2   3   4   5   6   7   8   9  10  11  1 2  13  14  15 


1   0   0 

0   1   1 

IN  PAIRS,  THE  NUMBER  OF  K-EOGEO  PATHS  AND  BRANCH  LOAOING 

EDGE  12   3  4   5  "~6  1      8~  9  10  11  12  13  14  15 


I   1   0 
1   1   1 
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RUN  NO.      1   FROM     1   TO 


Q-SUB-I,J  =  C. 9809998 


TffF VALUE  OF  Q-SUF-TT,J)  WITHOUT  BRANCH   1 

q-sub-i,j  =  o.aQ9-9aaaL — _ 


THE  VALUE  OF  Q-SUB-(I,J)  WITHOUT  BRANCH   ? 
Q-SUB-1,J  =  C.  8  9999  99" 


THE  VALUt  OF  C-SUB-dtJ)  WITHOUT  BRANCH   3 
Q-SUB-I ,J  =  0.8999999 


THE  VALUE  OF  U-SUB-- (  I  ,  J4  WI TH  BRANCH     I 


HAVING  RELIABILITY  VALUEO . 500000 
Q-St*8— I, J  =  0.9049999 
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THE    GIVEN^PATH    WAS 

2       3       4       5       6 

100     a — cr 


EOGE    1 

0 


7       8       9    10    11     12    13    14    15 
~0 — TT   ~0~    0       0       C       0      OF 


MATRIX  OPC  FOLLOWS 

~~r"   3   4   5" 


EDGE  1 

e- 

1 

C_ 

1 
1 


1 

0 

1 

0 

1 

1 

1 

1 

0 

1 

8   9  10  11  12  13  14  15 


THERE  ARf     2  PATHS 
THE..  PATH  MATRIX  IS 
EDGE  1   2   3   4   5   6 


8       9    10     11     12    13    14    15 


IN    PAIRSt     THE    NUMBER    OF    K-EDGED    PATHS    AND    BRANCH    LOADING 
2       3       4       5 — 6 T    8  ~9    10    11     12    13    14    15     ' 

4 1_&~  — _ _ . 
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RUN  NO.      2       (-POM     1   TH 


)-SUB-l,J  =  0.9309993 


THE    VALUE    OF    0- SUB- TT7JT ""WITHOUT    BFANCH    '1 

.  Q-suB-i.j  =  o.8S9iaaa 


THE  VALUF  OF  Q-SUB-(UJ)  WITHOUT  BRANCH   2 


Q-SUB-I ,J  =  C.R09999P 


THE  VALUE  OF  G-SUB-(I,J)  WITHOUT  BRANCH   3 
Q-SUB-I, J  =  0.  3999999 


THE- -VALUE  OF  Q-SUd^i-UJ)  WITH  dRANCH     2 
HAVING  RELIABILITY  VALUEO. 5OOOG0 

Q-SUB-I, J  *  ^.9049999 


MATRIX  QK  FOLLOWS 

1  1   C   — - 

1  0   1 

1  1   1 

0  1   1 


Q-SUB-K  =  0.9719998 
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THE  GIVEN  PATH  WAS 
EDGE  12   3   4   5   6   7 

0   0   1   0   C   0   0 


R   9  10  11  12  13  14  15 

0   C   0   0   C   0   0   0 


MATRIX  QPC  FOLLOWS 
FDGE  12   3   4   5 


c 

0 

1 

1 

0 

1 

n 

1 

1 

l 

1 

1 

l 

1 

0 

9     10  11  12  13  14  15 


THERE  ARF     2  PATHS 
THE  PATH  MATRIX  IS 
EDGE  12   3   4   5   6 
C   0   1 


10  11  12  13  14  15 


IN  pairs,  THF  DUMBER  OF  K-EDGEQ  PATHS  AND  BRANCH  LOADING 
EDGE  1   2   ?   4   5   6   7   8   9  10  1 1  1 2  1 3  14  15 


1   1 


&2 


RUN  NC.      ^       FROM     2   TO     3 

0-SUB-ItJ  =  C.98C9998 

THE  VALUE  OF  0-SUfw  ( I , J )  WITHOUT  BFANCH   1 

Q-SUB-I,J  =  C. 8 399999 

THE  VALUE  OF  Q-SUB-C I« J I  WITHOUT  BRANCH   2 

0-SUB-I  ,J  =  i.8^'9<?^ 

THE  VALUE  '.IF  G-SUB-UtJ)  WITHOUT  BRANCH   3 

Q-SUB-I ,J  -    C. 8099999 


THE  VALUE  OF  (,-SUti- (  I  » J  )  WITH  BRANCH     3 
HAVING  RELIABILITY  VALUEO. 500000 

Q-SUB-I, J  =  0. 90*9999 


MATRIX  OK  FOLLOWS 

1   0  1 

0  1  1 

1  1  1 
1   1  G 

O-SUB-K  =  0.9719998 
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THE  TOTAL  NUMBER  OF  PATHS  FOUND  IN  EACH  RUN 

—a—*—* —  - 


IN    PAIRS,     THE    NUMBfcK    OF    K-EDGED    PATHS    AND    BRANCH    LOADING 


RUN    NO.  1       FPQM  1       TO 

1    ~l~~"o~ 
1       1       1 


RUN    NO.  2       FROM  1       TO 


1       1       C 
1       1        1 


RUN    NO.  3       FROM  2       TO 

1       1       1 


TOTAL    BRANCH    LOADING 
3       3       3 
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